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ABSTRACT 


This  report  treats  the  theory  of  the  stability  of  an  elastic  rotor.  It  covers 
such  effects  as  instability  caused  by  friction  damping  in  the  rotor,  instability 
caused  by  shaft  asymmetry,  instability  caused  by  asymmetry  of  the  rotor  mass, 
instability  induced  by  fluid  film  Journal  bearings,  the  effect  of  static  damping 
on  the  stability  of  the  rotor,  the  influence  of  gyroscopic  effects,  the  effect 
of  flexible  bearing  supports ,  and  the  effect  of  gravitational  and  unbalance 
forces . 

The  report  gives  the  basic  governing  equations  and  the  methods  of  their  solution. 
The  conditions  are  established  under  which  instability  is  encountered  and  it  is 
shown  how  the  results  relate  to  a  practical  rotor  system. A 
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displacement*,  both  linear  and  angular,  fro*  the  uli  of  rotation.  Many 
of  ch<  result*  Are  well-known  end  era  summarized  hare  Ju*t  for  the  sake 
of  completeness ,  Some  result*  are  believed  to  be  new,  or  At  least  were 


new  when  first  obtained  by  the  author. 

The  restriction  to  email  displacements  renders  the  dynamical 
equations  of  motion  linear.  This  linearity  allows  one  to  use  super¬ 
position,  thus  computing  separstcly  the  rotor  motion  ' '.rider  -f facts  of 
gravity,  rotor  unbalance,  and  the  free  motion  of  the  rotor,  that  1$, 
its  motlou  near  equilibrium,  in  absence  of  any  external  forces  and 
torques.  Most  of  the  following  Is  confined  to  a  study  of  the  free 
motion.  The  affects  of  external  forces,  however,  era  considered  briefly 
In  Section  H. 


Am  regards  the  rotational  speed  of  the  rotor,  it  la  assumed  throughout 
most  of  the  following  that  this  spaed  Is  maintained  fpatfflt •  this  may 
require  a  special  drive,  such  as  an  electric  motor,  turbine,  or  belt 
drive.  In  order  to  simplify  mettors,  the  mechanism  for  insuring  the 
constancy  of  speed  Is  neglected,  as  is  also  the  frictional  torque. 

However,  the  effect  of  variable  rotor  spaed  Is  briefly  considered  In 
Section  11,  largely  In  running  through  the  critical. 

With  gravity,  unbalance  or  other  external  forces  absent,  the  differ¬ 
ential  equations  for  the  "free  motion"  of  the  rotor  become  both  linear  sad 
homogenous.  When  tha  coefficients  of  these  differential  aquations  are 
constant,  their  solutions  consist  of  linear  combinations  of  axponantlmls 


(1) 


where  the  *.’•  ere  roots  of  a  proper  algebraic  equation,  known  as  tha 
"characteristic  equation": 


\n  +  a^  i.11’1  + 


If  iU  roots  \  .  ..\a  o£  (2)  have  a  raal  negative  part,  than 


-1* 


«vtry  aolutlon  of  t.  dyna*U.cel  lyttM  representing  displacements  fro* 
•quilibriu*  approach**  zero  *•  time  inert****.  On  the  other  hand,  if 
env  one  of  the  root*  of  (2)  he*  e  positive  reel  pert,  then  the  general 
•olutlon  inert****  indefinitely  with  tint,  end  the  equilibrium  poeltlon 
1*  unatable. 

From  point  of  vice  of  (tabillty  it  i*  not  really  neceaeary  to  aolve 
the  characteristic  alfebraic  equation  (2)  for  all  of  it*  root*;  it  1* 
sufficient  to  determine  whether  each  root  has  a  negative  real  part: 

fce  (1^)  <  0;  1  B  1>  2,  ...»  n  ,  {3) 

In  tha  following,  such  equations  will  be  referred  to,  quite  inaccurately, 

as  "stable  equations".  Likewise,  we  refer  to  (2)  at  a  "real  aquation", 

when  its  coefficients  a.,  a.,  ...,  a  are  all  real  numbers,  or  can  be 

u  i  o 

made  real  by  dividing  (2)  by  a  constant.  An  aquation  which  Is  not  "real" 
in  the  above  sense,  will  be  referred  to  as  a  "complex  equation".  s 

A  well-known  test  for  the  stability  of  a  real  aquation  is  due  to 
Adolf  Hurwltz  (laf.l).  It  is  : 

Let 

eo>0.  (4) 

From  the  coefficient  of  (2)  form  the  n>fn  matrix 


(5) 


a 


2n-l  **•• 


a 


n 


where  the  subscripts  of  the  elements  in  the  first  column  increase  in  steps  of 


2  from  one  row  to  the  next  one,  while  the  subscript*  of  the  element*  in  any 
one  row  decresse  regulsrly;  when  the  subscript  is  lsrger  then  that  of  eny 
coefficient  occurring  in  equation  (2),  or  negative,  that  sleaent  is  re- 


*i  1  m  Us*  *«**»**• 
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*^  ■  0  for  k<0  or  k>  n 


(6) 


A  necessary  and  sufficient  condition  for  the  stability  of  the  real 
equation  (2)  is  that  the  following  inequalities  hoi  1: 

h>  °>  V0'  (7) 

all  hold,  where  ,  ..  0q  are  the  "principal"  determinants  of  the  matrix 
(5),  given  by 


■l* 

*1  *0 
*3  *2 

*1  “o 
a,  a 


0 


3  *2  *1 


*5  *3 
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D  -  a.  a  0 

n  1  o 


I *2n-l  •  *  ‘  1 

An  alternative  stability  criterion  for  a  real  equation,  dua  to 
South,  (Ref.  2)^111  also  be  recalled. 

Break  up  £(  \ )  into  two  parts  as  follows: 
f(\)-U(\)+V(>0 

wheve 


U(  K  )  ■  aQ  \n  +  *2  ^  +  .  . 

V( \)  -  a^0-1  +  \V'3  +  - 


(10) 


where  0  contains  All  the  tiru  of  (2)  with  even  eubeeripta  (that  ie,  the 
Ctiu  with  «o,  t2 ,  . ..)»  whlla  V  contains  the  terns  with  odd  subscript*. 

ZL-riim  7  lots  U,  denoting  toe  quotient  oy  q  enc  tM  remainder  by  Li 

C  (  N)  «  7U)  Q<  \)  +  K(  k).  <U) 

Then,  In  order  chat  the  real  equation  (2)  be  ■table  it  ia  oaceaaary  and 
aufflclent  that 

ao  ‘l  ?  °*  (12) 

and  (if  n>l),  the  equation  of  degree  (n  -1): 

fj  <  \)  ■  v<\)  no  (13) 

be  stable.  For  n»l  (12)  la  both  neceaaary  and  aufflclent. 

Jhtrthar  applications  of  the  above  lauth  criterion  enable  one  to  reduce 
the  condition  for  stability  to  a  number  of  inequalities  of  the  form  (12)  and  to 
the  stability  of  an  ever  lower  degree  equation. 

Turning  to  complex  equations  (2) ,  two  teats  for  stability  will  now 
be  given. 

▲  neceaaary  and  aufflclent  condi  for  the  stability  of  a  complex 
equation  (2)  ia  that  the  real  aquation  of  dagrea  2tu 

f(\>  f(\)  -  [ao\n  +  ..  +  ej  [io **•+...  +  in]-  0  (14) 

ba  atabla.  Bara  bars  denote  conjugates: 

c  +  id  ■  c-  id  {15) 

and  f  (\)  la  obtained  from  f (>■)  by  replacing  all  tha  coefficients  of  the 
lattar  by  their  conjugates. 

A  test  for  complex  aquations,  which  ia  similar  to  the  touch  teat  fives 
above  for  real  equations,  la  aa  follows. 

Let  eo  in  (2)  be  real,  end  denote  the  following  complex  coefficients  thee] 


•k  "  \  +  1  ck’  co  “  0 
Break  up  f (k)  a*  in  (9) ,  where  now 

D(\)  -  boXn  +  i  c^*'1  +  b2>.n‘2  +  t  c3\“'3.+, 

V(\)  -  bn\n‘l  +  t  c2\n'2  +  b3>.n"3  + . 


d«) 


(17) 


DividiV  into  U,  u  in  (11),  denote  tha  quotient  by  Q  end  the 
remainder  by  X.  men,  in  ordar  that  the  complex  equation  (2)  be 
•table,  it  is  neceeaary  and  aufflclent  that 

b  b.  <  0  ,  (18) 

0  1 


and  that  the  complex  equation  (13),  of  degree  n-1,  be  stable. 


The  leat  theorem  wee  proved  by  the  authot  la  an  Internal  General 
Electric  report, ’stability,  Hunting,  Follow-Ofc,  Issued  In  1938,  where 
various  other  testa  for  stability  were  also  presented.  After  a  lapse 
of  a  somewhat  long  interval  (26  years),  it  Is  planned  to  publish  soon 
the  rather  interesting,  purely  geometric  proof  of  It. 


In  certain  cases  tha  linearised  rotor  dynamical  equations  possess 
variable,  periodic  coefficients,  generally  of  a  period  equal  to  half  the 
period  of  rotation.  It  Is  known  that  special  solutions  of  these  equations 
exist,  of  the  form 


Pt(t) 


(19) 


vhare  are  periodic,  of  the  same  period  as  the  coefficients  and  are 
constant.  The  stability  of  the  solutions  now  reduces  to  the  requirement  (3) 
for  the  exponential  multipliers. 


Unfortunately  no  simple  test  similar  to  the  above  -riteris  for  etability 
of  the  solutions  (19), of  a  differential  equation  with  periodic  coefficients 
exists.  These  equations  will  be  briefly  discussed  in  Section  X  . 


5- 


Consider  a  rotor  coneleting  of  a  diak  of  aua  a,  eymeetrically 
■upported  on  a  relatively  light  ahaft  of  atiffneaa  k,  and  rotating 
with  conatcnt  angular  valoclty  about  a  flxad  axia  00  (taa  Flg.i), 

Aaauaa  no  uobalanca,  naglact  frictional  torquea,  or 

d  -  I/k 

llMiUX  1 

rathar  aaauaa  that  any  frictional  torque  la  ovarcoaw  by  a  proper 
driving  torque,  ao  that  a  conatant  epeed  of  rotation  la  aalntainad. 

Denote  by  a,  y  tha  displacement  of  tha  rotor  (and  abaft)  cantor , 
relative  to  fixed  (l.a. Newtonian)  (*,»,*)  -  axes,  with  origin  on  tha 
c-axla,  at  the  center  of  the  untleflected  ahaft.  Then  tha  equations 
of  the  rotor  notion  are 

a*  *  -kx  ,  my  <•  -ky.  (1) 

In  the  above  vs  have  neglected  forcea  other  than  the  one  due  to  the 
ahaft  atlffnesa  k. 

The  solution  of  (1)  leads  to  slaply  harmonic  motion  of  x,y,  of 
(radian)  frequency 

u>  -  ^k/o  (2) 

This  will  be  called  tha  "critical  speed"  of  the  rotor.  More  explicitly 

x  -  A  coa(oJot  +  «)  ,  y  -  B  cos^t  +  s')  (3) 

If  the  phase  angles  «,  «'  are  aqual,  then  the  path  la  a  straight 
line  through  the  origin.  For  general  e,  the  path  auy  be  ahown 
to  be  on  ellipse.  If 

c'  ■  e  +  n/2,  A  -  4  B  (4) 

the  peth  it  e  circle. 
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V*  n*xfc  add  to  the  abaft  fore#  a  '’frictional"  fores,  opposing 
tbs  velocity: 

»  -  - 

-  fit  ,  -  fy  ,  (5) 

where  f  la  a  positive  constant.  This  fores  vs  taro,  somewhat  In¬ 
appropriately,  "static"  friction,  in  tha  sense  that  it  opposes  the 
velocity  relative  to  static,  l.e.  fixed,  axes.  The  equations  of 
translational  notion  of  the  rotor  center  now  bacons 


nx  «  -  kx  -  fx, 
mJC  ■  -  ky  -  fy. 

The  solution  of  each  one  consists  superposition  of  exponentials 


(6) 


where  \  satisfies  the  equation 
m\*  +  £\  +  k  ■  0. 

Tha  two  roots  of  (8)  are  given  by 

- 


\ 


4  ka 


2  a 


(7) 


(8) 


(9) 


and  each  one  has  a  negative  real  part.  Thus  the  displacements  fron 
the  undeflected  position  x  ■  0,  y  ■  0  approach  aero  with  tine;  the 
affect  of  friction  la  to  dampen  out  any  initial  disturbance  fron  the 
equilibrium  position,  which  position  is  thus  seen  to  be  a  truly  stable 
position. 


In  the  above  it  will  be  noted  that  the  rotation  of  tha  rotor  has  no 
effect  on  the  bending  of  the  shaft  and  the  notion  of  the  rotor  canter. 

In  fact,  Bqs.  (1),  (6)  and  (9)  do  not.  involve  co,  and  they  remain  the 
same  whether  tha  rotor  is  rotating  or  not.  The  rotation  of  the  rotor 
is  merely  superposed  on  the  top  of  the  displacaaent  of  the  rotor  center, 
the  shaft  rotating  about  its  bent  center  line,  the  rotor  and  shaft  center 
undergoing  the  displacement  x,  y  corresponding  to  a  force  equal  to  the 
right-hand  member  of  (6)  applied  at  its  canter. 
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In  general,  S‘  i*  email  compared  with  the  t *rm  ‘m,  and  the  roots 
K  in  (9)  may  be  put  in  the  fora 


2a 


±  i 


-\  1 ^  ■ 


<10) 


where  i  •  - 1 , it  the  reel  pert  of  X,  end  ±  the  imaginary  part, 
leading  to  expressions  for  both  x  end  y  as  linear  combinations  of  products 
of 


-v 


±i\jt 


(11) 


The  pure  imaginary  exponentials  represented  by  the  second  factor  in 
(11)  when  properly  combined,  yield  real  sinusoidal  tame  of  frequency  X 


given  by 


J 


f I  -  (h>2  -  ■  [l  ‘  2  4»>  ‘  8  4/  +"]  <l2> 


while  the  real  exponential  in  (11)^  as  stated  above,  leads  to  an  ex¬ 
ponential  !ecay  of  amplitude.  For  f  -  0, i.e.  in  absence  of  friction, 
reduces  to  the  critical  speed  <dq  given  by  (2),  while  for  small  f. 


"j  ’“o(l  •  Is ). 


(13) 


The  introduction  of  friction  thus  leads  to  a  slight,  second  order 
decrease  in  the  frequency  Nevertheless,  we  shall  continue  to  refer 

to  (oq  (rather  than  X.j)  as  the  critical  speed  of  the  rotor. 

For  f  ■  0  the  motion  has  been  described  above.  For  f  >  0,  the  general 
motion  may  be  viewed  as  describing  an  ellipse  whose  axes  remain  fixed 
in  direction,  but  whose  semi-axes  are  decreasing  exponentially.  For 
special  motions  the  ellipses  become  circles  and  then  the  path  is  a  logar  ithmic 
spiral;  these  are  represented  by  the  equations 

'■  -\2fc  iM 

x+iy-Ae  ze  j  v 


+  iy  -  B  e"X2te“iXjt * 


(14) 

(15) 
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The  general  motion  say  also  be  visvsd  es  §  superposition  of  the 
two  logartbaie  spiral  motions »  given  by  £qs.  (14)  end  (15);  else 
as  obtainable  from  either  (14)  or  (15)  by  effecting  a  homogeneous 
strain,  or  akewing^of  the  (x,y)  plane,  of  the  type  in  which  psrallal 
straight  linea  go  into  parallel  straight  lines,  bus  distances  end 
angles  ere  not  preserved. 


Summarizing,  it  has  been  shown  above  that  the  effect  of  static 
friction  on  e  balanced  rotor  consist;  iu  reducing  eny  initial  dis¬ 
placement  of  the  rotor  center  to  zero,  both  the  x  and  the  y  coor¬ 
dinates  varying  with  time  as  exponentially  damped  sinusoids.  The 
motion  of  the  rotor  center  takes  place  in  an  ellipse  whose  dimensions 
decrease  exponentially  with  time,  the  frequency  being  given  by  X.^, 
the  damping  being  determined  by  given  by  equation  (10).  Generally 
is  close  to  the  critical  frequency  cdq,  and  \r  la  independent  of 
the  rotational  speed  of  the  rotor.  , 


Ill 


•  -  uswt-m&mQ.  m„  PESttfc  jeomi 

The  elastic  shaft  force,  for  a  round  shaftV  is  directed  along 

a  radius ,  from  x,  y  toward  the  center  of  the  undeflected  shaft;  it 

may  thus  be  described  as  a  "radial"  elastic  (return)  force.  To 

this  elastic  force  Y  in  2  (1)  we  now  add  a  force  F  ,  normal  to  the 
'vr  ^,n 

displacement  (x,y)  and  given  by 

in  ^'V  *  klx)  (1) 

where  k^  is  a  constant  (see  Fig. 2). 


Figure  2 


The  equations  of  motion  2  (1)  are  replaced  by 

e  • 

mx  ■  -kx  -  k.y, 

(2) 

my  ■  -ky  +  kjX. 


In  the  system  (2)  the  variables  x,  y  do  not  separate  as  in  2  (1). 
This  system  can  be  solved  by  putting 


x  •  X  e 


\t 


y  ■  Y  e 


\t 


(3) 


leading  to  an  algebraic  equation  in  of  degree  4.  A  slight  simplifi¬ 
cation  results  by  proceeding  as  follows: 

Multiply  the  second  Eq.  (2)  by  i  ■  FT  and  add  to  the  first  one.  In¬ 
troducing  the  variable 

Z  ■  x  +  1  y 


(4) 


on*  1*  ltd  Co  eht  atcond  ordtr  differential  equation 

aZ  -  (-k  +  1  kt)Z  (5) 

Thi»  It  aolvtd  by  putting 

Z  -  A  aXt  (6) 

resulting 

m\^  +  (k  -  1  k^)  ■  0  , 


2 

Since  A.  complex,  one  o£  it*  two  root*,  aay  will  have  a  real 

positive  part,  the  other,  a  real  negative  part  (see  Fig. 3). 


The  general  solution  of  (5)  is 

Z.Al6¥+A2e-V  (8) 

where  A^^  are  arbitrary  complex  constants.  The  A^-term  (i.e.  the  term 
involving  A^)  becomes  infinite  with  time;  the  A^-term  approaches  zero. 

From  (8)  x  and  y  are  obtained  by  taking  real  and  imaginary  parts  respectively. 

Thus  the  effect  of  the  force  F  normal  to  the  displacement  and  pro- 

portional  to  it,  is  to  make  the  rotor  motion  unstable.  This  conclusion 
holds  whether  k^  is  positive,  as  in  Figs.  2,  3,  or  negative,  that  is, 
whether  the  normal  force  is  in  the  direction  of  rotation  or  in  the 
opposite  direction. 


The  path  described  by  the  A^  -tens  on  the  right  of  (8)  re¬ 
presents  motion  in  a  logarithmic  spiral,  with  |zj  increasing  ex¬ 
ponentially  with  time;  the  second  term  represents  motion  in  a  log¬ 
arithmic  spiral,  of  the  same  angular  frequency  but  rotating  in  the 
opposite  direction  and  j  z\  decreasing  exponentially  with  time.  The 
general  solution  (8)  is  the  vector  resultant  of  the  two  motions,  and 
as  time  gets  large,  approaches  the  receding  logarithmic  spiral. 

The  notion  of  a  normal  force  is  not  purely  an  academic  one.  For 
a  steam  turbine  or  a  gas  turbine,  bending  of  the  rotor  shaft  de¬ 
creases  the  radial  clearance  between  the  rotor  and  stator  blades  on 
the  side  of  the  displacement.  This  causes  a  decrease  in  the  gas  leak¬ 
age  over  the  edges  of  the  rotor  and  stator  blades,  with  a  consequent 
Increase  in  the  local  gas  reaction  on  the  rotor  blades.  On  the  opposite 
side  an  opposite  effect  takes  place.  However,  the  Increase  in  force 
on  the  displacement  side  generally  overbalances  the  decrease  on  the 
opposite  aide.  Thus,  there  results  a  net  (resultant)  force  F  ,  normal 

S*  ti 

to  the  displacement  of  the  center. 


It  will  be  shown  in  Section  IX,  that  the  lubricant  in  a  journal 
bearing,  generally  also  produces  a  force  F  on  the  Journal,  normal  to 
the  displacement  (in  the  bearing)  of  the  Journal  center. 


Returning  to  Eq(8)  it  is  will  be  noted  that  the  unstable  term  Aje 
eventually  leads  to  increasing  energy  of  rotor  motion.  This  energy 
can  only  grow  from  the  work  W  done  by  the  force  F_  +  F  ,  during  the  rotor 
motion,  where 


If  k1/k 


kx  -  k^y)  x  +  (-ky 

- T +  + 

i  J. 


+  ^  x)y]  dt 

yx)  dt 


is  small,  then 


Icd 

o 


0) 

o 


(9) 


(10) 
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V 

and  the  exponential  growth  of  e  per  cycle  is  small, and  so  Is  the 
decay  of  e  ^ t .  Then,  for  the  e^lC  term  in  (8),  the  integral  in 
(9)  per  cycle  becomes 

kLr2  /  d^tan"lj)-  2x  k^2  -  2fl  kjxj  2  (11) 

Por  the  e^lC  it  becomes 

-  2  k2s  |a2  |  ;  (12) 

while  for  the  cross  terms  the  integrals  vanish. 

It  is  clear  from  the  above  that  if  |a2|  in  ( 8  )  la  much  larger 

than  |  AjJ  ,  then  the  rotor  motion  will  be  largely  clockwise  and  the 

force  Fn  will  take  energy  out  of  the  rotor.  Eventually  the  Aj-term 

will  become  equal  to  and  exceed  the  A, -term  and  then  F  will  feed 

2  ^.n 

energy  into  the  motion. 

If,  on  the  other  hahd,  the  direction  of  F  were  reversed  in  Fig. 4 

**  X  t 

by  replacing  k.  by  -k..then  F  would  cause  decay  of  the  e  1  -term 

1  -Art  '-n 
and  growth  of  the  e  1  -term. 

This  explains  why  the  motion  is  unstable  irrespective  of  the 
sign  of  k^.  Each  sign  promotes  its  own  exponential  motion. 
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IV 

ROTARY  FRICTION  FORCE 


We  now  return  to  the  ceae  considered  in  Sectionll,  but  replace 
the  static  friction"  force  2 (5)  by  a  "rotary  friction"  force.  By 
this  i a  meant  a  force  proportional  to,  but  opposite  in  direction  of, 
the  velocity  of  the  rotor  center  relative  to  a  rotating  system  of 
axes,  rotating  about  the  (fixed)  z-axis  with  the  same  (constant) 
angular  velocity  to  as  the  rotor.  This  relative  velocity  may  be  ob¬ 
tained  by  subtracting  from  the  actual  velocity  (x,y)  the  entrainment 
velocity  due  to  the  rotation  of  the  system,  given  by  (-ay,  ax)  (see 
Fig.  4).  Denoting  the  coefficient  of  rotary  friction  by  f  ,  we  now 
obtain  for  the  equations  of  motion  replacing  2(6). 

mx  -  -kx  -  f^(x  +  ay), 

(1) 

my  -  -ky  -  f1(y  -  ax), 
where  a  is  the  velocity  of  the  rotor. 


Introducing  Z  as  in  3(4),  one  obtains 

mZ  +  fjZ  +  (k  -  ifxa)Z  -  0.  (2) 

Once  more  this-  can  be  solved  as  in  3(6)  where  X  is  now  the  root 
of  the  complex  equation 
2 

mX  +  fjX  +  (k  -  ifjOQ)  -  0. 

To  insure  stability  each  root  X  must  have  a  real  negative  part. 


(3) 


'VXWWBM***** 


For  a  3econd  degree  complex  algebraic  equation 

b0X2  +  <bj  +  lcx)\  +  (b2  +  ic2)  -  0  (4) 

where  b  ,  b.  ,  .  cn  are  real,  the  teat  deacrlbed  in  the 

text  accompanying  Eqs  .  1(16),  1(17)  leads  to  the  following  In¬ 
equalities,  \ 


bo  bl  >  °’ 
<blb2  +  clc2) 


boc2 


>  0, 


(5) 


as  a  necessary  and  sufficient  condition  for  the  stability  of  (4). 


Applying  (5)  to  (3)  one  is  led  to 

mfl  >  0,  (6) 

,f.  (k  -  mco^)  >  0. 

w  *y 

The  first  inequality  in  (6)  is  automatically  satisfied  (since  f^ 
is  positive).  The  second  one  leads  to 


co  <  k/m  ■  “o'  (7) 

Thus,  for  rotor  speeds  below  the  critical,  the  rotor  motion  is  stable. 


For  rotor  speeds  greater  than  the  critical  - 

CD  >  CO 

O 


(8) 


the  second  inequality  in  (6)  will  fail,  at  least  one  of  the  roots  of 

\ 

(3),  say  will  have  a  real  positive  part,  and  the  rotor  motion  will 
be  unstable.  Since  the  sum  of  the  two  roots  is 


X1  + 


(9) 


S; 

V 


! 

i 

i 

! 

i 


i>'8 
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it  follows  that  for 


R*(\)  >  0,  **(X2)  <  0 


The  transition  to  instability  occurs  for 

03  ■  £0  a  V  — 

o  J  m 

when  X^  is  pure  imaginary.  Eq.  (3)  now  yields 

m  Xx2  +  lc  -  0, 

f.X,  -  if, CD 
11  1  o  , 

and  hence 


X,  «  Icd  • 
1  o 


(10) 


(ID 


(12) 


(13) 


Thus  the  transition  to  instability  occurs  when  the  rotor  is  running 
at  the  critical  speed,  and  it  is  characterised  by  a  circular  synchronous 

whirl  in  the  direction  of  rotation,  with  the  shaft  "frosen"  in  its 
whirling,  deflected  shape. 

It  f ollow8  from  (9)  that  the  root  Xj  l®3^®  to  an  oppositely  rotating 
whirl,  which  dampens  out  with  time. 

1 

In  the  unstable  range  (8),  the  root  X^  yields  an  exponentially  in¬ 
creasing  whirl.  While  the  unstable  range  is  independent  of  f^,  it  can 
be  shown  that  the  rate  of  growth  of  unstability  increases  with  f^. 


The  notion  of  "rotary"  friction  is^ itself,  by  no  means  academic. 

Rotary  friction  is  produced  by  any  device,  such  as  rubber  pa<h  or  rubbing 
surfaces  mounted  on  the  rotor.  Shrink  fits,  especially  of  wide  disks 
or  plates  shrunk  on  the  shaft,  are  another  cause  of  rotary  friction,  since 
if  the  shaft  bends  the  shrunk  element  does  not  follow  the  axial  extension 
and  compression  due  to  bending  of  the  shaft.  Thus,  slippage  takes  place 
along  the  shrunk  surfaces  when  the  shaft  bending  changes^and  the  shaft 
center  displaces  relative  to  a  system  of  axes  rotating  about  the  e-axis 
with  the  rotor  speed. 


Likewise  "solid  friction"  of  the  sheft  material  contributes  to  the 
rotary  friction;  by  "solid  friction"  is  implied  the  small  hystersis 
loop  obtained  by  plotting  the  strain  vs  stress  for  a  periodic  stress- 
strain:  even  in  the.  elastic  range,  the  curve  is  not  a  straight  line 
but  forms  a  thin  loop  whose  area  is  a  measure  of  the  energy  lots  per 
unit  volume,  per  cycle.  For  stresses  within  the  elastic  limit,  the 
energy  loss  per  cycle  is  of  the  form 

f  V  (14) 

where  V  is  the  maximum  elastic  strain  energy  and  f  is  a  "solid  friction" 
coefficient,  characteristic  of  the  shaft  material.  For  larger  stresses 
f  is  not  constant,  but  increases  with  amplitude . 


It  is  not  evident  that  the  above  rotary  friction  effects  are  truly 
represented  by  the  f^-terms  in  (1),  namely  as  a  "viscous"  force,  pro¬ 
portional  to,  but  opposing  the  relative  velocity.  The  following  consid¬ 
erations  show  in  fact  that  they  are  not.  „ 

For  sinusoidal  strains  and  displacements,  solid  friction  effects 
are  more  correctly  treated  by  replacing  the  Young's  modulus  E  by  the  com¬ 
plex  number  / 


(1  -  i  s)  E, 


i-pr 


Thus,  for  the  complex  solution  3(4),  if  time  is  assumed  to  enter  as  a 
factor  . 


Eq.  3(2)  can  be  treated  by  replacing  the  shaft  stiffness  constant  by 

k(l  -  i  s)  (17) 

and  this  is  analogous  to  Eq.  2  (6  )  if  the  friction  terms  -  fx,  -£y  are 
replaced  by  -kipx,  -  kipy  corresponding  to  a  value  of  s  in  (15)  pro¬ 
portional  to  the  frequency: 

8  »  kp  .  (18) 

It  is  of  interest  to  note  from  Eqs.  (1),  (2)  that  the  "rotary  friction" 

force  terms  involving  f^  can  be  represented  as  the  sum  of 

-  f  Z  +  ifx  mZ  (19)  . 

of  a  static  damping  force  and  a  force  perpendicular  to  the  displacement, 
and  analogous  to  Fn  in  3(1)  but  with  a  constant  k^  proportional  to  the  rotor 
speed:  k  „  f  m  (20) 
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Vs  now  consider  a  rotor  powaftaaiag  both  static  and  rotary  friction  in  a 
manna r  similar  to  that  of  Sections  II  &  IV.  The  aquations  of  motion  era 

not  +  fx  +  f  j  (x  +  ay)  +  lex  »  0, 

my  +  fy  +  f£y  -  arc)  +  ky  ■  0, 

Introducing  Z  as  in  3(4) ,  ona  obtains 

mS  +'  <f  +  f,)  i  +  (h  -  ifx  od)  Z  «  0. 

One  is  led  to  solutions  3(6)  where  \  are  roots  of 
m  \2  +  (f  +  ft)  \  +  (k  -  if j<o)  •  0. 

Applying  the  criteria  4(3)  for  stability  Of  this  complex  equation,  one 
obtains  the  inequalities! 

m(f  +  f.)  >  0, 


(f  +  fj)2*  -  mfx2  oj2>0. 


The  first  inequality  is  again  automatically  satisfied,  while  the  second  one 
can  be  put  In  the  dimensionless  form 

(l+f-)2  >  <S-)2-.(B-)2 

loo 

leading  to  !'-v 


|->a-  -l. 

£1  “o 


Two  cases  arise.  If  the  rotational  speed  is  less  than  the  critical  speed  <dq, 
t'-en  the  right-hand  member  of  (6)  is  negative,  and  (6)  is  automatically 
satisfied  irrespective  of  what  the  values  of  f,  f^  are.  On  the  other  hand, 
if  cd  is  greater  than  cdo,  the  right-hand  member  of  (6)  is  positive,  and  (6)  will 
be  satisfied  provided  that  the  static  friction  coefficient  f  is  sufficiently 
large  compared  with  the  rotary  coefficient  f^. 

The  relation  between  cd/<d  and  f/f,  is  shown  in  Fig.  5  where  the  solid  line 

o  i 

AB  marks  the  boundary  between  the  region  of  stability  and  instability.  We  thus 
arrive  at  the  criterion  that  if  the  rotor  is  running  above  its  critical  speed, 
there  will  be  stability  provided  that  the  static  friction  be  sufficiently  large 
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compared  to  the  rotary  friction  or,  what  amounts  to  the  earn  thing,  that  the 
rotary  friction  be  sufficiently  snail  compared  to  the  static  friction. 

As  an  example,  for  ao/ao^  ■  2  corresponding  to  a  rotor  speed  which  is  double 
the  critical  speed,  f/f^  must  exceed  unity,  that  is,  the  coefficient  of  static 
friction  f  must  at  least  be  equal  to  the  coefficient  of  rotary  friction  f^,  to 
insure  stability.  For  to/o^  -  2,  if  f^  >  f,  there  will  be  instability; 
if  ^  <  f,  stability  will  obtain. 

In  the  unstable  range  the  rate  of  increase  of  whirl  increases  with  fj/f. 

Finally,  if  to  the  forces  in  (1)  one  also  adds  the  normal  force  3(1),  one 
is  led  to 

mx  ■  -kx  -  k.y  -  fx  -  f . (x  +  coy), 

(7) 

my  ■  -ky  +  k^x  -  fy  -  f^  (y  -cox). 

Introducing  Z  as  in  III  (4) ,  one  is  led  to 

raZ  +  (f  +  fj)  Z  +  (k  -  i  k^  -  i  ©)  Z  -  0  (8) 

and  to  exponential  solutions  for  Z  as  in  3(6)  where 

m\  2  +  (f  +  fL)\  +  k  -  i(kt  +  oa^)  -  0.  (9) 
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By  applying  similar  techniques,  one  now  obtain*  for  •tablllty  of  (9)  the 
following  inequality 

f,  „„  k-fl) 

U  >  a  .  l  +  _l_2  do) 

£i  “o  h 

The  boundary  of  the  region  of  stability  i*  shown  a*  the  dotted  slanting  line  A1  B* 
on  Pig.  5,  for  positive  kjj  for  negative  k^  A'B'  ia  to  the  right  of  AB. 

Remarks  similar  tc.  those  made  at  the  ends  of  Secs  III,  IV  apply  equally 
well  here.  In  some  cases  it  is  possible  to  fit  the  above  theory  so  as  to 
Include  effects  of  solid  friction  and  other  aspects  of  frictional  losses,  by 
assuming  that  the  "constants"  k^,  f,  f^,  vary  with  cu  (or  even  with  the 
amplitude  of  oscillation).  When  this  is  done,  it  is  evident  that  the 
rectilinear  boundaries  of  the  stability  region  of  Fig.  5  may  be  greatly 
distorted. 
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VI 

p^mrnu  mrr  .mrmss  mmms. 

We  now  consider  a  roLor  aounted  on  a  shaft  of  uniform  but  non-circular 
cross-section.  A  normal  section  of  such  a  shaft  possesses  two  mutually  per¬ 
pendicular  directions  corresponding  to  the  maximum  and  minimum  moments  of 
inertia  and  it  possesses  different  bending  stiffness  constants  k^,  kg*  kg'*  ^1’ 
in  these  two  mutually  perpendicular  directions.  One  may  define  two  vibrations 
frequencies 


As  a  result  of  the  two  stiffnesses  it  can  be  shown  that  a  force  in  a 
direction  different  from  the  principal  directions  of  inertia  produces  shaft 
deflections  lying  in  a  plane  different  from  the  plane  through  the  z-axie- con¬ 
taining  the  force.  To  find  the  displacement  one  resolves  the  force  along  the 
principal  directions  of  the  shaft  section,  determines  the  deflection  along 
each  one  and  superposes  them. 

If  the  shaft  is  rotating  and  the  force  is  fixed,  then  the  deflection  com¬ 
ponents  turn  out  to  involve  in  the  coefficients  sln2tot,  cos2tot,  where  o>  is 
the  velocity  of  rotation  of  the  shaft. 

As  a  result  of  the  above  the  differential  equation  of  motion  also  involves 
the  above  trigonometric  functions. 

To  avoid  this  complication  we  refer  the  motion  to  a  system  of  axes,  (x,y) 
which  rotates  about  the  fixed  z-axis  with  the  rotor  (and  shaft)  speed  to.  The 
deflection  of  the  rotor  center  is  given  by 

xi  +  yi  (2) 

where  i,  j  are  unit  .vectors  along  the  (rotating)  x,y-axes. 

These  i,j  are  not  constant,  but 

AJ 


di 

dj 

-X  =  cn/i  ’ 
dt  J 

d?  "  €  “jb 

(3) 
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Differentiating  (2)  we  have  for  the  (true)  velocity 

di  di 

v-il  +  yj  +  x^f+y^:,  W 

end  substituting  from  (3) 

v  «  (x  -  coy)  i  +  (y  +  cox)  j  (5) 


Further  differentiation  yielda  for  the  acceleration 

a  m  0c  -  2coy  -  to2*)i  +  +  2cdx  -  <e2y)  j.  <S) 

Hence  the  dynamical  equations  of  the  rotor  are 


*  2 

m(*  -  2fl)y  -  cd  x)  -  -kj  x  * 

•  2 

m(y  +  2 cd ^  x.  -  cd  y)  -  k2  y . 


(7) 


This  system  has  constant  coefficients.  The  variables  x,  y  in  (7) 
do  not  separate.  Hor  can  these  equations  be  simplified  by  intro¬ 
ducing  complex  variables  as  in  3(4). 


Assuming  a  solution  of  the  form  3(3)  one  is  now  led  to  the  fourth 
degree  equation 

(\2  -  CD2  +  CDj2)  (\2  r.  CD2  +,  C^2)  +  4<D *\2  -  0.  (8) 

2  2 

This  is  quadratic  in  \  .  Hence,  unless  both  roots  \  are  negative, 
the  solution  .  Eq.J  (7),  )»111’  ho  unstable. 


Put  (8)  in  the  form 


where 


yf*  +  2x3>+  c  ■  o, 

(9) 

b  -  (CDX2  +  cd22)  /2  +  cd2>0 

2  2  2  2 
c  -  (cd^  -  CD  )  (cd2  -  CD  J* 

.(10) 

The  \  roots  of  (9)  are 

\2  -  -  b  ±\/(b2  -  c) 
The  radicand  above  is 


b2  -  c 


9  2.2 

“l2  .  “2  \  2  2  2. 

L  ~  +  7ix>  (cd.  ^  +  c or)  >  0, 


(ID 

(12) 
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an*  i*  always  positive.  Hence  the  roots  (11)  are  always  real.  If 
c  is  positive,  then  both  roots  (11)  are  negative.  On  the  other  hand, 
if  c  is  negative,  then  the  root  -b  +{7~-c  will -be  positive,  and  8C53S 
solutions  of  (7)  will  be  unstable.  Examination  of  (10)  shows  that 
this  happens  for 

oq12  <  oj2  <  <o22  (13) 

Thus,  for  a  rotor  with  a  shaft  possessing  two  different  stiffness 
constants  in  two  mutually  perpendicular  directions,  the  speed  range 
(13)  between  the  two  vibration  frequencies  (1)  is  unstable. 

If  static  friction  also  exists,  then  the  range  of  instability 
narrows,  and  for  a  sufficiently  large  value  of  the  coefficient  f, 
may  disappear  altogether.  On  the  other  hand,  normal  forces  and 
rotary  friction  may  enlarge  it. 

Thus  far  we  have  considered  only  translational  displacements  of 
the  rotor.  Angular  displacements  of  the  rotor  axis,  it  will  be  shown 
in  Section  VIII,  introduce  gyroscopic  effects,  and  may  produce  a  further 
unstable  range  at  high  enough  speedb. 
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VII 


GYROSCOPIC  EFFECTS 

Even  when  the  rotor  is  mounted  as  in  Pig.  1,  it  may  execute  both 
translational  notion,  normal  to  the  a -axis,  ae  veil  es  rotational  notions  about 
transverse  axes.  The  latter  motions  are  even  more  likely  to  occur  whan  the 
rotor  is  mounted  asymmetrically,  for  Instance  if  it  la  overhung,  as  indicated 
schematically  in  Fig.  6.  The  following  analysis  applies  to  anv  mode  of 
support  of  the  rotor. 


x 


Let  the  c-axls  contain  the  undeflected  shaft  center  line  poaltlon, 
end  denote  by  x,y  the  normal  translational  displacements  of  the  rotor  canter, 
end  by  \,  p.  its  angular  (or  rotational)  displacements  about  the  x,  y-axes; 
ell  four  quantities  x,  y,^,n  are  assumed  to  be  smell. 

In  this  Section  we  assume  that  the  shaft  is  axially  symmetric  elastically, 
and  that  tha  rotor  is  axially  symmetric;  in  Section  VIII  ve  consider  the  case  vhen 
the  shaft  stiffness  constants  in  two  mutually  perpendicular  directions  ere 
different,  and  the  rotor  has  different  moments  of  inartls.I,  J,  X  4  J, 
about  the  principal  exes  of  inertia,  normal  to  the  spin  axis;  the  third  prin¬ 
cipal  moment  of  Inertia,  about  the  spin  axis,  will  be  denoted  by  K. 

For  the  axially  symmetric  shaft  (and  rotor),  a  force  F^,  applied  to  the 
shaft  at  the  rotor  center,  and  a  moment  (about  the  y-axis),  produce  a  linear 
displacement  x  and  a  rotation  n  about  the  y-axis,  given  by 


x  -  a  F  +  b  T  , 

x  y 

U  -  b  F  +  c  T  , 

x  y 


<D 


where  a,b  are  the  linear  and  angular  displacements  produced  by  a  unit  force  F 
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acting  alone;  b ,c  tha  similar  displacements  produced  by  a  unit  torque 

about  the  y-axie ,  acting  alone.  The  matrix  of  coefficients  a ,b,c  in  (1)  is  the 

■’flexibility  matrix".  Solution  of  (1)  yielde 

Fx-  Ax  +  B4.  (2) 

T  ■  Bx  +  CH, 

where  the  matrix  A,B,C  ia  the  "etiffnees  matrix",  end  is  the  inverse  of  the 
matrix  in  (1). 

Similarly,  an  applied  force  end  moment  T^  (about  tha  x-axis)  produce 
displacements  in  the  (y,  z  ) -plane,  gir'en  by 


y  - 

X  ■ 


a  F 


-b  F 


b  T 


c  T 


(3) 


and  by  the  inverse  equations 

Fy  *  Ay  -  B  X  ,  ^4) 

-  -By  +  C  X  . 

If  the  e-axis  is  a  stable  position  of  equilibrium,  then  the  quadratic  form 
Q  -  (Ax2  +  2  Bx  U  +  C  H2)  /2  (5) 

v 

representing  the  elastic  energy  corresponding  to  the  deflection  (x,  is 
positive  definite.  The  conditions  for  this  are 

A>0,  AC-B2>0,  (C>0)  (6) 

(where  C>0  follows  from  the  preceding  conditions).  Similar  remarks  apply  to 
the  matrices  (1),  (3),  (4). 

It  is  assumed  above  that  the  displacements,  both  linear  and  angular,  ere 
small,  so  that  superposition  applies,  and  the  order  of  the  rotations  X,  n  Is 
immaterial. 

It  is  further  assumed  that  the  rotor  is  balanced  about  its  axis  of 
rotation;  this  assumption  implies  that 

(a)  the  center  of  mass  (c.m.)  of  the  rotor  lies  on  its  axis  of  rotation; 

(b)  the  axis  of  rotation  is  a  principal  axis  of  inertia  of  the  rotor, 
through  its  c.m. 

These  conditions  are  assumed  to  obtain.  The  assumed  rotational  symmetry 
implies  equality  of  the  two  transverse  moments  of  inertia,  I  »  J.  Hence  any 
two  axes  through  the  c.m.  mutually  perpendicular  and  normal  to  the  s-axis,  can 
be  considered  as  principal  axes  of  inertie  of  the  rotor  in  its  undeflected 
position.  In  particular  the  unit  vectors  i,  k,  parallel  to  the  x,  y,  x-axes. 
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can  be  used  as  principal  directions  of  inertia,  even  though  they  are  fixed 
in  apace  and  hence  not  in  the  rotor.  After  the  (small)  angular  deflections 
the  above  principal  directions  vectors  become: 


il  -  i 
jl  " 


- 

J  +i^, 

•v  v 

i|i  -f  +  k  . 


(7) 


Noe  let  vary  with  time.  The  motion  of  the  rotor  R,  at  any  time  instant, 

consists  of  a  translational  velocity  of  its  c.m. 

v  -  ii  +  yj  ^  (8) 

^  -v-» 

and  of  a  rotation  represented  by  a  proper  rotation  vector  through  its  c.m. 
given  by 

cd  -  \i  +  |lj  +  cok.  (9) 

assuming  that  R  continues  to  rotate  about  its  bent  shaft  with  velocity  cd. 

Utilising  Eq.(7),  cd  can  also  be  put  in  the  form 

•  •  v 

cd  -  1.1 .  +  p.)  +  cbk.  (10) 

provided  second  order  terms  in  ji:  jl  are  neglected.  This  form  yields 

the  component  of  cd  along  the  (Instantaneous)  principal  axes,  and  this  simpli¬ 
fies  the  expression  of  M,  the  moment  of  the  momentum,  to  be  referred  as 


m.m.: 


M  -I  (Xij^  +  HJ^)  +  KtDkj 

a  ^ 

-I(\i  +  jij)  +  Kcd(1ji  -  J\  +  k) 


(ID 


where  the  last  form  is  obtained  again  by  using  (7)  and  neglecting  second 
order  terms . 

Differentiation  of  M  yields 


M 


(1  \  +  Kn  CD)  i  +  (I  H  -  K\  CD)  j . 


(12) 


We  now  proceed  to  set  up  the  dynamical  equations  of  motion  of  R.  Recall 

that  in  Eq.  (1)  F  i  -4-  F  j  is  the  force  applied  to  the  shaft  at  the  c.m.  of  R; 
* 

this  force  is  either  an  external  force  or  it  is  applied  by  R  to  the  shaft. 


■ttWH 


The  reaction  force,  -F  i  -  F  J,  la  therefore  applied  by  the  shaft  to  R. 

x~  y  ~j 


Similarly  the  shaft  applies  the  torque  -T^i  -  T^j  to  R 


Equating  these 

shaft  reactions  on  R  to  the  mv  and  H,  we  obtain  the  dynamical  eqvstions 


mx  +•  Ax  +  Bji  “0, 

•• 

my  +  Ay  -B\  -  0, 

•  •  • 

I\  +  Kyui  -By  +CX  -0, 

—  9 

1^  -  KXtD  +Bx  +Cn  ■  0. 

Assuming  a  solution  of  (13)  in  which  time  enters  as  a  factor 
vt 


(13) 


(U) 


;nd  rearranging  the  order  of  the  equations  there  results  the  following 
determinantal  equation 

_ X _ M _ 2 _ 


A(v)« 


A  +  m  v 
B 
0 
0 


B 

e 

C  +  I  vi 
0 

K  ojv 


0 

0 

A  +  m  v" 
-B 


0 

-K  co  v 
-B 

C  +  I  v' 


(15) 


-  o. 


where  the  coefficients  in  any  one  column  arise  from  the  variable  indicated  in 
that  column,  in  the  row  above  the  determinant. 

The  function  A(v)  is  evidently  even  in  v,  and  of  the  fourth  degree  in 

2  2 

v  .  Thus,  unless  all  four  roots  v  are  negative,  at  least  one  root  will  have 
a  real  positive  part  and  the  solution  of  (13)  will  be  unstable. 

We  expand  A(v)  in  a  Laplace  expansion  in  terms  of  its  first  two  columns: 

A  -  (12)  [34]  -  (13)  [24]  +  (14)  |23] 

+  (23)  [14]  -  (24)  [I3j  +  (34)  (12}  (16) 

where  (ij)  is  the  2nd  order  determinant  in  the  first  two  columns  and  rows  i,j 
while  [ijj  is  the  determinant  in  the  last  two  columns  and  the  same  rows. 

There  results 


A  -  (1-2)2  +  (14) 2 


(A+mv2)  (C+Iv2)  -B2J  2+^(A+mv2)  K  V  (a  j  2  «  0 


(17) 


ry 


Bsssusaas® 


(18) 


To  consider  negative  v2,  we  put 


.  2  2 
ip  *  v  -  -p  . 


Then  A  can  be  factored  thus 

A  - 


(A-mp2)  (  C-Ip2)  -  B2]-  (A-mp2)  K  p  m 


(A-mp2)  (C-Ip2  )  -  B2]+  (A-mp2)  K  p  <o 


}• 


(19) 


The  first  factor  is  linear  in  a>  and  when  equated  to  sero  yields 

2  - 

2>  1 


CD  -  f(p) 


1  2.  B2 

'  K "  l  '*Vl 


(20) 


The  right  hand  member  is  odd  in  p,  and  has  the  asymptotes 

P  -  0,  p  -  1  jAM  ,  cd  -  -(I/K)  p  .  (21) 


Near  p  -  0 

cu  ■ 


(CA  ~  BZ)  1 
KA  p 


(22) 


where  is  analytic  at  p  ■  0.  If  we  assume  that  the  s-axls  is  a  position 
of  stable  equilibrium,  then jnoting  Eq.  (6)j  the  residue  at  p  ■  0  is  positive. 
Likewise  the  residue  at  p1  -"^A/m  is  positive. 


A  schematic  plot  of  Eq.  (20)  is  shown  in  Fig.  7-  It  i 3  evident  that 

a  line  of  constant  <o>0,  cuts  this  curve  in  4  distinct  points,, 

p  "  pl*  p2*  p3*  P4J>  <23) 

two  positive,  two  negative;  but  these  roots  cannot  be  negatives  of  each 
other.  They  correspond  to  the  time  factors. 


»lPfc 


P  "  Pj»  J  -  1*2»3»4- 


The  second  factor  in  (19)  leads  to 
o>  -  -f(p) 


(24) 


(25) 


and  hence  on  Fig.  7  to  the  reflection  of  the  curve  shown,  in  the  p-*axis, 
and  to  4  roots  which  are  the  negatives  of  (23). 
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It  has  thus  been  shown  that  all  four  roots  v  of  A  are  indeed  nega¬ 
tive  and  distinct.  Thus  for  symmetric  rotor  and  shaft  the  gyroscopic 
effects  never  lead  to  instability. 

The  factorization  of  A  suggests  that  the  system  (13)  may  be  reduces xe 
_o  two  systems  of  lower  order. 

By  multiplying  the  second  Eq.  (13)  by  i  and  adding  to  the  first  one,  the 
third  one  by  -i  and  adding  to  the  fourth  one,  one  is  led  to  the  second  order 
system  of  differential  equations 


mZ+AZ  +  BM  =>0, 

•  *  • 

B  Z  +  I  M  -  i  KxH  +  C  M  =  0, 


(26) 


where 


Z  =  7  +  iy,  M  -=  p  -  iX,  (27) 

A  conjugate  set  of  equations  can  be  derived  for  the  variables 


Z  =  x  -  i  y,  M  =  u  +  i\. 


(28) 


Eqs.  (26)  and  their  conjugate  equations  could  have  been  used  in  place  of 
Eqs .  (13). 

Indeed,  assuming  a  solution  of  (26)  in  which  time  enters  as  a  factor 
ipt  ... 

c--  yields  the  frequency  equation 


(A  -  rnp ) 


3  >,(.  +  KLup  -  ip  ) 

Leads  to  Eq.  (20). 


=  0 


Bach  of  tha  four  roots  of  (23)  result*  la.  *  solution  of:  th*.  i«?s 

*  +  iy  •  Z0  #ip£f  n  IX  «  K0  eipt.  O'?) 

Taking  real  part*  ona  may  Interpret  tha  saotioa  as  a  whirl  of  tha  rotsr. 

In  which  the  shaft  executes  a  rotation  of  frequency  p,  with  each  paint 
on  the  axis  describing  a  circle  and  the  exit  ee  a  whole  a  hyperboloid  of 
rotation.  This  motion  may  be  described  as  a  circular  whirl.  In  parti* 
cular  the  point  on. the  rotor  axis  at  a  distance  s  from  tha  origin  describes 
the  circular  path 

x  +  iy  +  (*  -  *0)  (u  -  IX)  -  |zo+  (e  -  sQ)  KQ]  eipfc  (31) 

where  z  «  zq  corresponds  to  the  rotor  c 

The  first  factor  in  (19),  put  in  tha  form 

(A  -  xnp2)  [c  -(I  -  K  ®)  p2J  -  B2  -  0,  (32)  . 

shows  that  for  aach  circular  whirl  Sq.  (32)  is  tha  same  as  the  equation  of 
vibration  of  a  non-rotating  rotor  whoaa  transverse  moment  of  inertia  I 
is  replaced  by 

i'-  I  -  K  <o>/p)  (33) 

In  certain  cases  this  may  lead  to  negative  I',  but  always  to  real 
frequencies  of  vibration. 

Similarly,  from  the  conjugate  equations  for  Z,  M,  one  obtains  four 
circular  whirls  corresponding  to  the  negatives  of  the  roots  of  (29). 

As  a  special  case,  suppose  that  the  shaft  is  rigid,  and  is  restrained 
toward  z-axis  by  springs  of  (angular)  stiffness  C.  Then  the  equations  of 
motion  corresponding  to  the  moment  of  momentum  equations  about  the  origin 
suffice  for  describing  the  dynamical  behavior,  and  they  lead  to 

IX  +  Kuio  +  CX  ■  0,  (34) 

14  -  kxd  +  Cm.  =  o. 


where 


now  I  is  the  moment  of  inertia  of  the  rotor  about  the  x-  and  y-axes. 


(35) 


One  is  led  Co  the  frequency  equation 

(I  p2  -  C)2  -  K  p2  co2  -  0, 

Solving  for  o>,  there  results 

m  ■  +  (I  p2  -  C)  /Kp.  (35) 

A  plot  of  u>  vs,p>0  is  shown  in  Fig.  8.  Again,  for  each  constant  co,  there  ass 
now  two  real,  distinct  values  of  p,  corresponding  to  two  circular  ^recessions 
of  whirls,  both.  positive!  there  .are  also  two  negative  one^  (»ot  shown). 


It  is  of  interest  to  note  that  if  the  z-axis  is  not  a  stable  equilibrium 

position,  then  the  residue  in  (22)  is  negative,  and  some  of  the  real  roots 

2 

of  (20)  may  be  lost,  leading  to  positive  V  and  unstable  solutions,  esp¬ 
ecially  for  small  u>. 

For  the  special  case  implied  in  Eqs.  (34)  such  an  unstable  equilibrium 
along  the  z-axis  corresponds  to  negative  C.  For  negative  C.  Fig.  8  is 
modified  into  Fig.  9,  with  a  minimum  for  co  at  the  point 


It  will  be  noted  that  for 

co;>cd  (38) 


that  are  two  real  root*  of  (36)  corresponding  to  two  precessions!  fre 
quancia*  and  stable  solutions  of  (34).  On  the  other  hand,  for 

a)  <<a  .  (39) 

o  > 

there  exist  no  real  roots  of  (36);  hence,  the  roots  of  (36)  are  com¬ 
plex  and  the  motion  Is  unstable. 


From  the  above  it  is  evident  that  gyroscopic  effects  may  stabilize  an 
otherwise  unstable  position  of  equilibrium. 

An  illustration  of  the  above  occurs  for  a  top  spinning  on  a  horizontal 
table.  Here  the  negative  C  arises  from  the  gravity  force.  For  a  high  enough 
speed,  the  top  "goes  to  sleep",  with  its  center  of  mass  above  its  point  of 
support.  As  the  rate  of  spin  decreases  as  a  result  of  friction,  the 

ft  n 

position  becomes  unstable,  the  top  wobbles,  and  eventually  falls  down. 


VIII 

STABILITY  OF  A  RQTATTOSALEY  ASYj-Sgl^SC  &OSSS  MOCSBSB' C5L.A  SSag?.KIS5.^ • 

DIFFERENT  STIFFNESSES 

Let  I,  J,  K  be  the  three  principal  moments  of.  inertia  of  the  rotor 
through  its  c.m,  K  being  the  accent  of  inertia  about  the  rotation  axis,  jfe  con- 
sider  a  system  S  of  (fc,  y,  z)  -  axes,  with  unit  "sectors  i,  j,  k,  rotating  re- 
lative  to  fixed  space  about  the  z-axis,  with  the  constant  rotor  velocity  <o. 
Suppose  that  in  the  undeflected  position  the  directions  i,  j correspond  to  the 

#*w  es/ 

principal  axes  of  inertia  of  the  rotor,  as  well  as  to  the  principal  axes  of 
bending  stiffness  of  the  shaft  crosssection,  that  is,  to  the  largest  and 
smallest  moments  of  inertia  of  that  section  about  axes  through  its  center. 

Suppose  that  relative  to  axes  of  the  rotating  system  S,  Eqs.  7(1)  -  7(4) 
hold  for  the  deflections  of  the  shaft,  but  different  sets  of  stiffness  constants: 


AjX  +  B^p,  Ty  -  B^x  +  C^p, 

(1) 

A2y  "  B2^’  Tx  ”  “B2y  +  C2^  * 

(2) 

Here  x,y;  \ ,|i  represent  the  (small)  linear  and  angular  deflections  of  the  rotor 

shaft  at  the  rotor  c.m.,  z  »  zq,  relative  to  the  rotating  axes  of  S,  while 

F  ,  F  ;  T  ,  T  represent  the  forces  parallel  to  these  axes  and  the  moments  about 
x  y  x  y 

these  axes  acting  on  the  shaft  at  the  c.m. 

Since  the  system  is  rotating  about  the  z-axis,  it  follows  that 


- 


0) 


l-  & 


~(x>  i. 


(3) 


To  avoid  ambiguity,  we  suppose  that  the  angular  displacements  X.  ,  p.  are 
carried  out  as  follows.  First  the  rotation  \  is  carried  out  about  the  vector  i; 
this  rotates  the  vectors  i,  j,  k  into  i.j  j. ,  k  ,  given  by 

i,  °  i,  j.  ■  j  cos  \  +  k  sin\,  k.  ■  -  j  sin\+  k  cos  \  .  (4) 

Next  carry  out  the  rotation  p  about  the  vector  leading  to 

i  =  i,  cos  n  -  k  sin  |i,  j0  »  j.  ,  k  =  i.  sin  p  +  k  cos  p.  (5) 

^  rS-  v  TV- 

Eqs.  (4),  (5)  apply  at  any  time  t,  with  the  values  of  p  corresponding  to  that 

time  instant.  The  vectors  i  ,  j  ,  k  correspond  to  the  princioal  directions  of  the 

1  rJ-  ^S- 

rotor  at  the  time  t,  in  its  deflected  (and  rotating)  position;  they  are  fixed  in 
the  rotor  R, 

From  Eqs.  (3),  (4),  (5)  follows  that  the  net  rotation  vector  of  R,  relative 

is  given  by 


(il  K  4  L,  x 


P-. 


(6) 


For  small  u  Eqs. 


(4) ,  (5)  reduce  to 


i,  - 1.  jt  - 1  +*»•  . 

h-h  -y-  ij-Ji-  i2-ii4+  -1- 

A<  ^ 

provided  pover.  ofh,  p  hlghar  then  the  fleet,  be  ne.lactad;  heoco 

i,  -  1  -  |ik,  3,  "  3  +  '"k,  kj  »  lit  -  *■  3  +  k. 

A«  ^ 

The  Inverse  of  Eqa.  (4),  (5)  C.ke  on  the  sene  form,  but  .1th  sin),  , 
.l„  ,  (or  4.  u)  repleced  by  their  nagatlvaa.  *>.  (6)  ^ 

of  term a  linear,  in  1,  only» 

*  • 

(D  -  cd  (k9  -  v  i,  +  * J2*  +  i2  X  + 

•  i  (X  -  <d  p)  +  J2(£  +  )  +  a3k2* 

,  „  r  u„  d  ik  leads  to  their  time  rates: 

The  cross  product  of  c»  by  i2,  J2»  K2 


<n 


di2 


d^2 


dt  "  °^2  ‘  (oA  +  *^2  ’  dt  -V 

"  ^  -  ^2  * 

Alternatively  equation.  (11)  can  be  obtained  by  differentiating  (9) 
and  utilizing  (3).  From  Eq.  (10)  folio.,  for  the  m.m.  of  * 

H.Il(^mn)+J  j20i  +  ®  +  Kcok2» 

-v  ^ 

and  substituting  from  (9) 

[IW  (K  -  I)  u<o]l+  ['>*«•«"  A]_)+Kaik. 


Differentiating  and  utilizing  (3),  there  results 

M.  *  f  l"K+  (K.  -  I  -  J)  H  ®  +  CK'J) 

v  5 


+  [j  £  -  (K  -  I  '  J)\  “+  (K  -  I)  “2  ^ 


(9) 


(10) 


d^2 

-<ci2+  (X  -  a*i)k2  ,  -  (H  +  ^H2  (11) 


(12) 


(13) 


(14) 


From  Eqs.  (1),  to 
-J?i  are  obtained: 


(14)  and  6(6)  the  following  dynamical  equations  of  rotor 


*  (V  -  2oy  -  ca  x)  -!■  AjX  +  "  0, 

•  *  *  2 

a  (y  +  2osc  *oy)  +  A^y  -  >.  K  G» 


I  \  +  (X  -  I  -  J)  pm  +  (K  -  J)  m2  \  -  B2y .+  C2  \»  0, 

J  *p*  -  (K  -  I  -  J)  Xm  +  (K  -  I)  ©2  p  +  BjX  +  Cx  p  «  0. 

Assuming  •  ■oluClon  of  Eqs.  (15)  of  the  fora  7(14)  ona  is  led  to  tho 
deterainental  equation 

P  y  *■ 


A(V)- 


V2  -  CD2)  + 

B, 


2  m  CD  v 
0 


B, 


-  2bkdv 


fj  v2  +  (K  -  I)  CD2  +  cj  0  (J  +  I  *  K)  CD  V 


(o(v2-cd2)  +  -  B2 


(5,5) 


(K  -  I  -  J)  CDV 


-  B, 


4v2  +  (K-J)  CD2  +  C^J 


(16) 


Expanding  as  in  7(16)  shows  that  A (v)  is  even  in  V  as  well  as  in  cd.  Thus 

2 

only  if  all  four  roots  v  are  negative,  will  the  solution  of  (15)  be  stable,. 

As  a  special  case,  consider 

cd  *  0,  (17) 

that  is,  the  case  of  a  non-rotating  rotor.  The  (x,y)-axes  are  now  fixed. 

Eqs.  (15)  separate  into  one  pair  of  (x,p)  equations  for  the  (static)  vibrations 
of  the  rotor  in  the  (x,z) -plane,  and  a  pair  of  (y,  )  equations  for  similar 

vibrations  in  the  (y,z) -plane.  The  resulting  deterainental  Eq.  (16)  factors 
into  two  second  order  determinants,  agreeing  with  Eq.  7(15)  for  cd  «=  0.  The 
four  roots  of  (16),  provided  the  z-axis  is  a  position  of  stable  equilibrium,  are 
negative  and  distinct.  By  continuity  follows  that  Eq.  (16)  is  stable  for  a 
proper  range  of  small  speeds: 

cu  <  cd'  .  (18) 

As  a  more  significant  check,  consider  the  special  case 

I  -  J,  A  »  A2»A,  -  B2  ■  E,  ci  "  c2  “  0-  (19) 

of  an  axially  symmetric  rotor  and  shaft  --  the  case  considered  in  Sec.  7. 

By  proper  manipulation  of  Eqs.  (15),  and  by  introducing  the  variables 

Z=*x+iy,K=p  ~i\, 

or;,:  obtains  the  fourth  order  system 


(20) 


**  *  2 

a(Z+2io>Z-  a  Z)  +•  A  2  *  BH  »  S5 

BZ  +  CM+IM  -(1C  -2  I)  i  o>  M  +  {£  -  1}  £D2  M  »  0. 

Assuming  time  to  enter  as  the  factor 

1  q  t 
e 

one  is  led  to  the  equation 

2 

A  -  m  (q  +  cd)  B 

2  2  "  0 
B  C  +  K  (o>  +  a  q)  -  X  (q  +  a>) 

It  will  be  noted  that  this  agrees  with  7(29)  provided  that  one  replaces 

CD  by  0),  p  by  q  +  <D, 

and  notes  that  factoring  out  e  *  ^  from  e  *  P  in  7(31): 

ipt  icot  iqt 
e  r  ■  e  e  n 


(21) 


(22) 


(23) 


(24) 

(25) 


corresponds  to  passage  from  fixed  to  rotating  axes. 

The  plot  of  q.  vs  cd  can  be  obtained  from  Fig.  7  by  means  of  the  transformation 

CD  ■  CD,  q  ■  p  -  CD,  (24) 

that  is  by  moving  the  points  of  the  curve  of  Fig.  7  at  any  constant  height  cd 
to  the  left  a  constant  distance  equal  to  cd.  Again  one  is  led  to  four  values 
of  q  for  each  cd. 

The  transfconnation  (26)  at  first  sight  appears  to  be  at  odds  with  the 
conclusion  that  A  (v)  in  Eq.  (16)  is  even  in  V.  However,  it  will  be  recalled 
that  the  fourth  order  systems(21),  7(26)  hold  only  for  the  variable  Z,  M.  Ti;? 
conjugate  systems,  for  Z,  M,  lead  not  to  (26),  but  to 

0)  ■  0),  q  *  p  +  cd,  (27) 

and  this  corresponds  to  shifting  the  reflection  of  the  curve  of  Fig.  7  in  the 
horizontal  axis>  a  distance  cd  to  the  right.  Geometrically^ different  portions  of 
the  curve  cd  =  +  f(p),  cd  >  0,  are  skewedj  cd  ■  f(p)  to  the  left,  cd  *  -f(p)  to  the 
right,  so  that  equal  and  opposite  v-roots,  still  go  into  equal  and  opposite  roots. 

Since  the  special  case  (19)  is  thus  always  stable,  it  follows,  by  continuity, 
that  the  system  (15)  is  stable  for  sufficiently  slight  rotor  and  shaft  asymmetries. 

We  turn  next  to  another  special  case,  namely  when  all  the  elastic  constants 
vanish: 


A ,  =  A  £  =  0,  B1  =  B2  =  0,  Cx  =  C2  =  0. 

leads  to  an  unrestrained  or  free  rotor. 


(28) 
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Sqs.  (15)  now  ssparat*  Into  (sr,y)-  nqiintton*  i*ich  l*ad  to 


<V2  -  oj2)  +4  <d2  V2  •  (V2  +  ffl2) 
and  to  (\,  li)-  aquation*  front  which  follows 


m  0 


<29> 


I  J  v4  +  jl  J  +  (K  -  I)  (K  -  J)]  a2  v 


+  (K  -  I)  (K  -  J)  «  -  0 


(30) 


It  1*  of  interest  to  point  out  that  Eqs.  (29)*  (30)  also  follow  fro®  (16) 

when  the  elastic  constants  do  not  vanish,  if  to  is  allowed  to  get  large,  provided 

2  2  2 
it  is  assumed  that  the  roots  v  of  A  (v)  also  get  large,  so  that  v  /to  remain 

finite.  This  assumption- is,  indeed,  valid  for  (29),  (30).  Hence  the  roots  of  (16) 

v“'  2  2  2 

behave  like  the  roots  of  (29),  (30)  for  large  to  .  In  the  (to  ,  v  )  -plane  the 

2  2  2  2 
ratios  v  /to  correspond  to  slopes  of  asymptates  of  the  curve  A(v  ,  CD  )  ■  0. 


The  repeated  roots 


2  2 
v  ■  -  to  , 


CO 

A  i  to  t, 


(31) 


te  fc,  which  represent  uniform  rectilinear 


lead  to  solutions  e 

motion  of  the  c.m.^when  described  in  terms  of  x,  y  of  the  rotating  system  S 
Eq.  (30)  can  be  put  in  the  form 


a  -  (K-I)  (K-J)/  IJ, 


x2  +  (1  +  a)  x  +  a  -  0,  x  -  ^2  , 
and  factors  into 

(x  +  1)  (x  +  a)  -  0 

thus  yielding 

v2/a2  -  -1  ,  v2/cd2-  -(K  -  I) (K-J)  /  IJ  -  -a. 


(32) 


(33) 


(34) 


It  will  be  noted  that  if  K  is  the  largest  or  the  smallest  moment,  of  inertia,  the 
roots  (34)  lead  to  stable  solutions;  on  the  other  hand  if 

I  <  K  <  J  or  J  <K<I  (35) 

then  one  root  of  (30)  leads  to  instability. 

The  above  results  for  (34)  agree  with  the  classical  Poinsot  motion,  where 

the  hodograph  paths  near  the  principal  axes  with  largest  and  smallest  moments  of 

inertia,  are  closed  curves  (on  the  elliproid  of  inertia),  but  form  two  self- 

2  2 

intersecting  curves  through  the  intermediate  axis  of  inertia.  (The  root  V  53  -u> 
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Ee  turning  to  general  values  of  the  elastic  constants ,  the  icccnbiiity  In 
ctw  (33)  i*pli«s  tfcet  th*  solution  of  (13)  la  unstable  for  large  anengh  to. 


K  la  maximum  or  minimum,  no  coneluaioo  can  be  drawn  regerding  th*  stability  of. 

2 

(15) f  except  that  for  large  enough  a  the  roots  are  all  stable. 

2  2 

It  la  of  Interest  to  find  the  values  Of  <d  for  which  v  *  0  1*  c  root  of  (16) , 

2 

•lnce  passage  fro*  stability  to  instability  may  occur  by  v  passing  free  a  negative 
to  a  positive  value,  hence  through  v  ■  0. 

If  V  »  0,  then  x,  become  constants,  and  eqs.  (15)  separate  into  an 

(x,  u)  -  pair  leading  to 

[-  »  a>2  +  )  [(K  -  I)  a>2  +  C1]  -  Bj2  -  0,  (3b) 

and  a  (y,  \  )  -  pair,  resulting  in 

[-  a  <d2  +  A2  ]  [(K  -  J)  cu2  +  cj  -  822  “  °*  (37) 

A  study  of  the  roots  of  (36),  (37),  under  the  assumption  that  Sht  z-axia-is  a 

stable  position  of  (static)  squilibriua,  lssds  to  the.  -On.;  o.U5i.vtV  .  %' r  'v Ji) 

2 

holds,  three  roots  cn  of  (16)  for  v  «  0  are  positive,  one  ncjar  ivc.  Vtien  (35) 
does  not  hold,  then  ell  the  four  roots  of  (16)  are  positive  for  K>  I,  t>  J,  or 
two  are  positive  and  two  are  negative.  Plots  of  A  (to,  v)  for  v  «  0  are  indicated 
schematically  in  Pig.  10. 


2 

Eq..  (16),  when  v.._7anded  tn  powers  of  v  ,  yields 
2  8  6 

A(oj,v)=m  I  J  v  +  (  )v  +  ...  +  A(co,  0)  »  0. 

The  product  of  the  roots  of  (16)  is  thus  equal  to 


(38) 


x*.  x*  xi  m*  m 

1  2  3  4 


.  T  .. 

4*<a,  u; 

■2IJ 


Of) 


2  2 

Bane*  follows  chat  for  a>  corresponding  to  nsgatiTa.  £(«  >  0)  la  Fig.  1A  not 
all  four  roota  of  (15)  can  ba  negative .  Thane  rotor  epae4  intervals  thua  cor- 
r*  snood  to  uaatabla  solutions. 


Tmmmmgm  xVCm  •Civility  tu  instability  u y  occur  not  opiy  fiunngri  r  ptfiia|  froa 

negative  to  positive  valuaa-tbdssL  occur  at  tha  ©2-roote  of  (34),  (37)-  but  al»o 

2  2 

through  tuc  nagatlva  roota  ,  V2  moating  and  raeadlng  at  right  aaglaa  In  tha 
complex  plana.  This  happens  at  simultaneous  roota  of  tha  two  aquation* 


A(cd2,V2)  »  0,  3&(a>2,  V2)/^2  ■  0 


(40) 


In  cone  lua  ion  .for  an  uasyamstrlc  rotor  and  abaft,  thara  always  ax  1st  unstable 
rotor  spaad  rangas. 

Aa  an  example,  wo  raturn  to  the  rotor  of  Sactlon  6  with  tha  two  ahafe  atiffnaaaaa 
kl*  k2  (corr*,Pon<ling  Co  A^ ,  Aj)  and  the  centrally  supported  rotor,  but  assume 
chat  tha  lattar  has  two  different  transverse  moments  of  lnartla  I,  J,  I  ^  J.  Xa 
terms  of  tha  constants  in  Kqs.  (1),  (2)  wa  have 

Bi  "  B2  "  0  (41) 

Bqs.  (15)  now  separata  into  x,  y  equations  whose  solution  has  bean  found  in 
Section  6  gnd  K,  u  equations  leading  to  the  determinant  aquation 


f(V2  a2)  -[jv2  +  (K  -  1)  co2  +  Cl]  [iv2  +  (K  -  J)  co2  +  cj 

2  2  2 

+  (K  -  I  -  J)  V  co  -  0 


(42) 


2  2 

Since  £  Is  quadratic  in  V  and  co  ,  Eq.  (42)  when 

2  2 

yields  a  conic  section  C.  For  large  co  ,  V  ,  the 
of  (30);  thus  C  has  real  asymptotes  whose  slopes 
numbers  of  (34). 


2  2 

plottsd  in  ths  (co  ,  V  )  -plane, 
roots  of  Eq.  (42)  approach  those 
are  equal  to  tha  right-hand 


-39- 


Thus  far  It  has  bean  assumed  that  the  I,J  axes  of  Inertia  of  the  rotor  are 
parallel  to  the  principal  directions  of  bending  of  the  shaft.  Suppose  now 
that  this  is  no  longer  the  case,  and  that  while  Eqs.  (3)  -  (14)  apply,  the 
principal  directions  of  bending  of  the  shaft  make  an  angle  7  with  the  principal 
directions  i,  j  of  rotor  moments  of  inertia.  We  shall  indicate  the  modifications 
in  the  dynamical  equations  (15)  arising  from  this  misalignment  of  principal  axes. 
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quAtions  (1),  (2)  Imply  an  elastic  potential  energy  Q  of  the  ahaft  given  by 
for  y  -  0) 


2Q  -  (At  x2  +  2BX  xn  +  Cx  »2)  +  (Ajy2  -  2B2  y\  +  C^2)  (46) 

ith  0  <  7  <  «/2,  the  energy  Q  is  obtained  by  replacing  (x,y;  X.,n)  in  (46)  by 

x' ,y ' ,n ' )  where  x'..|i'  are  the  linear  and  angular  displacements  along  (and 
bout)  the  principal  shaft  stiffness  axes  x',y',  where 

x'  *«x  cos  7  +  y  sin  y,  y'  *  -x  sin  y  +  y  cos  y, 

\  cos  y  +  n  sin  y,  li’  »  sin  y.  +  p.  cos  y. 

’placing  x .  p  in  (46)  by  x*,  . p',  then  expressing  the  latter  in  terms 

£  x,  . .,  p  by  means  of  (47),  there  results  for  Q  the  quadrative  fora 

29  *1J  Ul  Uj'  AlJ  'A-11 

lerej  to  make  the  summational  notation  available^u^ ;  i«=  1,2, 3, 4  have  been  intro- 
iced  for  x,y;  A,p,  and  A^  are  shown  in  Table  I  and  where  the  further  notational 
bbreviatione 

c  ■  cos  7,  s  "  sin  y  (49) 

ive  been  made. 


TABLE  I 


2  2  2  2  2 
■’Al  c  +  A2  C  ,  Au  -B,  c  +  B2  S  ,  A 


13 


A22  "  C1  f2  +  c2  ®2'  *23 


33 


(A1  '  V  80 *  ti4  "  <V  V  8C 


(Bj  -  B2)  sc,  A24  »  (C2-  Cx)  sq 


(Ax  *2  +  A2  c2)  A34=  -  (B]_  82  +  B2  c) 


'44 


2  2 
a  +  c2  c  ) 


e  dynamical  equations  replacing  (15)  are  now  obtained  by  replacing  the  force  and 
rque  terms  in  (15)  (i.e.  the  terms  involving  ...  e^)  respectively  by 
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*iar  w 


(50) 


dx  "  Saj  "  *11  *  +  *12y  +  *13*  +  *14** 

iy  ■  lu2  "  *21'  x  +  *22y  +  *23*  +  *24 11 

*  ‘  *31  x  +  *32y  +  *33*  +  *34'* 

’  *  *41  x  +  *42y  +  *43*  +  *44“  ' 

The  resulting  equations  lead  to  a  corresponding  determlnantal  equation  replacing 
(16),  but  now  with  no  vanishing  eleaents. 
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IX 


EFFECT  OF  JOUBMAL  BEARINGS  ON  mLMHH 

Thus  far  each  bearing  has  been  assumed  to  act  as  a"fixed"  support  for  the 
rotor  shaft  or  journal,  so  that  while  the  rotor  shaft  rotates  freely  inside  the 
bearing,  the  point  on  the  axis  of  the  shaft  (or  journal)  at  the  bearing  center 
was  assumed  to  undergo  no  transverse  displacements.  Even  for  pre-cospressed 
ball  or  roller  bearings,  this  assumption  is  satisfied  only  approximately,  due  to 
the  compliance  of  the  bearings.  For  cylindrical  journals,  however,  where  there  is 
definite  radial  clearance  between  the  journal  and  bearing  surfaces,  this 
assumption  will  never  do,  in  view  of  the  added  freedom  of  motion  of  the  journal 
center  in  its  "circle  of  clearance". 

The  position,  in  the  middle  plane  of  the  bearing,  of  the  journal  center 
0^  relative  to  the  bearing  center  0,  for  a  given  load  on  the  rotor  (such  as  the 
rotor  weight) ,  automatically  adjusts  itself  so  that  both  bearings  "carry"  this 
load;  that  is  0^  moves  to  such  a  position  that  the  pressure  developed  by  the 
lubricant  in  the  bearing  has  a  net  resultant  acting  on  the  journal  which  is 
equal  and  opposite  to  the  applied  load.  Basically  this  pressure  distribution 
is  due  to  the  shaft  "dragging  in"  the  lubricant  into  the  convergent  portion  of 
the  film. 

It  has  been  known  for  many  years  that,  at  times  ,fhe  rotor 
vibrations  apparently  maybe  induced  by  the  lubricant  in  the  bearings.  This 
can  be  demonstrated  in  a  laboratory  model:  by  temporarily  reducing  the 
lubricant  supplied  to  the  bearing,  the  vibrations  diminish.  For  a  horizontal 
shaft  and  reasonably  high  loads  these  instabilities  occur  at  speeds  considerably 
higher  than  the  critical  speed;  for  lightly  loaded  bearings  or  for  a  vertical 
shaft  they  may  occur  even  at  low  speeds. 

To  include  the  effects  of  journal  bearings  on  rotor  stability,  it  is 
evident  that  not  only  is  knowledge  needed  of  the  steady  state  bearing 
forces  corresponding  to  a  fixed  position  0^  of  the  journal  center,  but  also 
of  the  transient  forces  arising  from  the  motion  of  0^  near  its  equilibrium 
position. 

For  simplicity,  we  consider  only  complete,  360°  ,  self-acting  (that  is, 
non-pressurized)  journal  bearings.  We  suppose  that,  as  in  Fig.  1,  the  rotor  is 
mounted  at  the  center  of  the  shaft,  and  that  the  shaft  is  supported  at  its  ends 
by  two,  exactly  similar,  perfectly  aligned  bearings.  We  assume,  furthermore,  that 
the  displacement  001  of  the  journal  center  from  the  bearing  center  is  the  same  at 
each  bearing,  at  all  times.  Then  the  same  force 


<V  »,> 


a> 


will  be  exerted  an  the  rotor  Journal  by  tho  lubricant  la  each  bearing;  bora, 

•a  iuCwd  lu  «ij«  12 ,  MovCvt  tua  i««l«  1  wuafvwiuw  v«  £***£«»  ■mo««ow  fwlMt 

vbaa  directed  toward  the  baarlot  cantor  0,  while  la  the  forca  coaaoaiat  at 
right  anglea  to  00^ . 


A  aoctloo  of  tha  bearing  (with  an  exaggerated  radial  clearance  C)  io  ah own 
in  Fig.  12.  Tha  Journal  cantar  0^  can  evidently  ha  locatad  only  within  tha 
"clearance  circle^'  with  cantar  at  tha  bearing  center  O^and  of  radlua  C.  Bence 


OOj  -c  C  0  <  a  <  1  ,  (2) 

where  c  la  Known  aa  tha  "eccentricity’’;  «  ■  0  corrasponde  to  coincident  O,  0^; 

e  ■  1  to  actual  contact  of  journal  and  bearing  aurfacca. 

Sat  up  axaa  aa  in  Fig.  12,  with  tha  y-axle  in  tha  direction  00,  of  tha 
Journal  center,  ateady-atate  dlaplacaaMnt,  and  let  8  be  aeaaured  Iron  tha 
y-axia.  The  gap  ia  given  by 

h  •  C  (1  -  a  coa  8);  (3) 


It  la  evidently  even  in  6 . 


Figure  12 

Alio  shown  on  Fig.  12  are  the  forca  conponenta  Fr>  acting  on  the  journal. 
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- 1 

u  mIL  u  the  "attitude,  eagle” 

?  •  tan'1  It  /Tr.  <4> 

that  is  the  angle  bttwin  the  applied  load^and  the  displacement  vector,  00^. 


The  fcrces  Fr  ,  Ffc  art  given  by 
L/2  2x 

F_  -  a  /  d«  /  p  coa  9  d8, 
f  -L/2  0 


L/2  2s 

F '  «  -a  /  da  J  p  a  in  8  dfi, 

C  -L/2  0 


(5) 


where  p  la  the  preaaura,  a  the  radlua  of  the  Journal,  L  the  axial  bearing  length, 
a  the  axial  diatance  along  the  bearing  aurfaca  from  lta  canter  line,  and  la 
poaltlva  whan  dlrectad  radially  Inward.  The  preaaura  aey  be  found  by  integrating 
the  Reynolds  equation,  which  for  an  incompressible  lubricant  la 

IT [  J55  <S^> +  h  (!#>]  <‘> 

where  n  -  vlacoaity  of  lubricant, 

U  - a oo  ■  linear  velocity  of  journal  due  to  lta  rotation^ 

W  ■  h  »  normal,  (l.a.  radial)  velocity  of  Journal  aurfaca,  poaltlva  whan 
directed  toward  0^. 

W  vanlahea  If  0^  la  stationary,  othervlae  It  la  given  by 

W  -  -  (uf  coe  8  +  ut  ain  8),  ■  C«,  u£  ■  Caa  (7),  ^ 

where  ia  the  radial,  the  tangential  componenta  of  velocity  of  0^,  and  a  la 
the  angle  formed  by  00^  with  a  fixed  direction  (on  Fig. 12  a  ■  0).  Hence 


W  ■  -d  coe  8  -  C«a  a^n 


(8) 


Substitution  from  (3),  (8)  into  (6)  yielda 


45- 


This  equation  has  to  be  solved  so  that  p  is  periodic  in  G  and  satisfies  the 
boundary  condition 

p  *  p^  at  z  «  +  L/2  (10) 

at  the  bearing  ends,  where  p#  is  aabient  pressure. 

Solutions  of  (9)  have  been  obtained  under  the  assumption  that  p  is  inde¬ 
pendent  of  zr 


P  “  p(0) 


(ID 


an  assumption  that  is  valid  for  a  long  bearing  except  near  the  ends,  see 
references  (4),  (5).  Substituted  in  (5)  they  lead  to 


Fr  -  -  Ke  f^  ,  Ft  -  K(o)  -  2a)  f  t , 

provided  Fr  is  measured  positively  outward,  where 

3  2 

K  ■  12x  pa  L/C  , 


(2  +  €2)(1  -  €2)1/2 


(12) 

(13) 

(14) 


Consider  first  stationary  0^: 

k  *  0,  a  *  0  .  (15) 


Eqs.  (12)  yield 

Fr  =  0,  Ffc  *  Ito  ft(e). 


(16) 
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This  lap lies  that 

♦  -  x/2  (17) 

The  vanishing  of  Fy  for  a  stationary  0^  also  follows  from  the  sytmaetry  of  the 
film  thickness  about  0*0  end  0  ■  jc,  as  a  result,  the  pressure  Increase  over  the  con¬ 
vergent  half  of  the  film  is  matched  by  an  equal  pressure  decrease  over  the  divergent 
film  half.  In  equation  form: 

p(6)  +  p(-6)  -  2p(0)  -  2p(x);  (18) 

This  states  that  p(0)  -  p(0)  is  odd  about  0*0,  and  p (0)  -  p(s)  is  odd  about  8  **  n. 
(That  p(0)  ■  p(x)  follows  by  putting  6  ■  *  in  the  first  eq.  (18)  since  p(8)  is 
periodic  of  period  2x)  .  Indeed,  eq.  (18)  follows  from  eq.  (9)  for  4  *»  0,  by  noting 
that  the  right-hand  member  of  (9)  is  odd  in  6. 

Likewise,  the  coefficient  of  e  in  Eq.  (9)  is  even  in  0,  and  the  resulting  con¬ 
tribution  to  Fr>  Ffc  in  (5)  leads  to  a  vanishing  Ffc. 

The  same  considerations  of  evenness  and  oddness  of  p  carry  over  to  solutions 
of  Eq.  (9)  for  bearings  of  finite  axial  length  L.  Eqs.  (12)  are  valid  for  them,  too, 
but  with  a  different  K,  and  with  f^,  ffc  functions  of  e  and  L/2a. 

Actually,  the  above  prediction  (17)  does  not  agree  with  observed  positions  of  0^. 
For  a  vertical  load  of  different  magnitudes  and  at  different  speeds  cd,  the  equilibrium 
positions  of  in  its  circle  of  clearance,  tend  to  cluster  near  a  locus  schematically 
shown  in  Figure  13,  and  approximating  a  semi-ellipse.  At  most,  (17)  holds  only  for 
small  e,  that  is  for  light  loads,  or  very  high  speeds. 


Figure  13 

The  explanation  of  the  above  discrepancy  between  theory  and  tests  is  due  to  the 
negative  pressures  which  are  allowed  in  the  solvation  of  Eq.  (9),  but  which  cannot 
be  realized  in  a  liquid  lubricant.  If  the  pressure  is  decreased  in  oil,  it  will  tend 
to  cavitate  and  release  dissolved  air  and  oil  vapor.  Therefore,  Eq.  (9)  must  be 
solved  subject  to  the  condition  ' 

P  >  °>  (19) 

and  putting 

p  =  0  over  (20) 
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where  A  denote#  the  cavitating  area.  Many  shape#  of  A  say  be  found, 
c 

A  'binimum  lo#8M  principle,  Cameron  and  Mood,  reference  (6  ),  have  proposed  the  con¬ 
dition  of  vanishing  normal  pressure  gradient  at  the  edge  of 


^  «  0  at  boundary  of  Ac ,  (21) 

On  Figure  12  the  cavitating  area  A  is  shown  schematically  between  the 
broken  lines  on  the  divergent  film  half.  With  (20)  holding  on  Ac»  Eq.  (18) 
fails  to  apply  for  at  m  0,  €  “0,  and  a  non-rero  F  results. 

For  a  long  bearing,  if  Eq.  (11)  is  assumed,  A£  is  bounded  by  two  lines 

of  constant  0,  and  these  are  determined  by  Eqs.  (20),  (21).  For  finite  L,  the 

determination  of  A  may  require  trial  and  error  methods  and  special  programing 
c 

on  high-speed  computing  machines. 

Physically,  the  boundary  A  may  vary  with  p^  in  Eq.  (10)  and  with  the 
nature  of  the  liquid  lubricant,  its  vapor  pressure,  and  its  dissolved  gases. 

While  (9)  remains  linear  in  p  outside  Ac>  the  introduction  of  Ac  whose 
boundary  may  change  with  the  right-hand  number  of  (9),  renders  the  problem 
non-linear,  and  superposition  of  solutions  no  longer  applies. 

It  will  be  noted  that  while  the  right-hand  member  of  Eq.  (9)  involves  cu, 
a  and  e,  the  two  former  occur  in  the  combination  to- 2ct-  Therefore,  the  integra¬ 
tion,  even  with  the  cavitating  conditions,  Eqs.  (19)  -  (21),  can  be  confined  to 

a  m  °> 

and  the  results  used  to  obtain  Fr>  Ffc  for  any  a,  by  replacing  as  with 

-  2a  m  ajg,  g  «  1  -  2a/aj.  (23) 

Assuming  a  to  vanish,  and  using  dimensional  analysis  considerations,  one 
obtains,  for  a  particular  ratio  of  L/2a, 

F  »  -ipcnf  (e.e'),  F  u'o  f  (e,e')f  (24) 

r  r  t  l 

whs-* re  f  f  are  dimensionless  measures  of  F  ,  F  ,  l  has  the  dimensions  of  a 
r  ’  t  r  c  . 

ler.rr'n  so  that-^u^has  the  units  of  a  force,  and  €  has  been  rendered  dimensionless 


by  introducing 


r  -  a* ,  € '  -  k/ai  *  daj/dT. 


It  any  be  shown  that  the  length  l  can  be  put  in  the  fora 


.  3  _  ,„2 

i  *  s  L/C  t 


as  in  Eqs.  (12),  (13).  It  follows  now  that  for  C£  t  0, 


Fr  *  ~lV-  (<u  "  2a)  fr  (e.^'/g), 


Ft  -in(<o  -  2a)  ffc  (e,  €  *  /g) • 


g  -  1  -  2a' 


For  small  c  cavitation  is  not  likely  to  occur.  If  L/a  is  large,  Eqs.  (12)  - 
(14)  may  form  a  fair  approximation.  Neglecting  powers  of  c,  e  higher  than  the 
first,  they  yield, 

Kcue 

Fr  “  -Ke,  Ffc  *  K(co  -  2a)  e/2*  -y  -  Kcite ,  (28) 

Note  on  Figure  14,  that  Ce,  Cea  are  the  radial  and  tangential  velocity  components 
of  0^.  Hence,  the  force  exerted  by  both  bearings  on  the  shaft,  in  rectangular 
components ,  is 

Fx  -  -2f£  -  fowj,  Fy  »  "2t|  +  ftug,  f  -  |  (29) 

where  are  the  components  of  00^  in  the  x-  and  y-  directions. 

Suppose  now  that  the  shaft  is  infinitely  rigid.  Then  the  equations  of  motion  are 


m  £  =  F  *  -f£  -  four),  my  ■  F  *  -2ff)  +  fui£ 

x  y 


Introduc ing 


C  -  I  +  in 


md  assuming  time  to  enter  in  the  form  e  '  ,  one  is  led  to 


(31) 


(32) 


mX2  +  2f K  +  ifco  -  0, 

an  equation  which  by  the  tests  of  Section  I  is  always  unstable. 

Suppose  now  that  a  radial  bearing  force  be  added  to  the  above  so  chat  (29)  is 
changed  into 

Fx  -  -  2fg  “ f con  -  Fy  -  -  2fij  +  f*$  -  kby  (33) 

where  kfe  is  the  radial  bearing  "spring  stiffness"  constant.  Equation  (32)  becomes 

2 

m\  +  2f\  +  (^  +  ifco)  »  0.  (34) 

This  equation  agrees  w_th  5(3)  provided ,  in  the  latter,  one  puts 


(35) 


(36) 

that  is,  for  speeds  below  double  the  critical,  defined  as  the  frequency  of  vibra¬ 
tion  of  the  rotor  mass  against  the  bearing  (radial)  stiffness.  For  all  higher 
rotor  speeds,  the  rotor  motion  is  unstable. 


f  S  f 


k 


V 


Hence,  applying  5(6),  we  find  that  the  rotor  is  stable  for  speeds 

a>  <  2(0*,  , 


ax 

o 
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Suppose  next  that  the  shaft  is  not  infinitely  rigid  but  possesses  a  finite  stiff¬ 
ness  constant,  k.  Then,  to  the  given  center  deflection  (£ ,  T|)  is  added  a  further 
deflection  of  the  rotor  center,  given  by 


x  «  Fx/k,  y  «*  Fy/k, 


(37) 


resulting  in  a  more  complex  set  of  dynamical  equations  than  (30) .  The  exponential 
solutions  e^*"  of  these  equations  are  determined  by  the  cubic 


3  ,  ,  b  +  k  ,  ,  „  2  .  k  .  ki 

2\  +  ( — - icu)  X  +  2\  —  +  — 

K  £  mj 


zk, 

K-f >  ' 


I  (%  * 


It  is  shown  in  Reference  (  7 )  that  this  equation  is  stable  for 


(38) 


co  <  2cuq 


(39) 
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and  unstable  for  higher  rotor  speeds,  where  a>Q  is  defined  by 


% 


“b 


<k  +  kj* 


(k  +  Kjj)  m' 


(40) 


This  critical  is  the  vibration  frequency  of  the  rotor  mass  against  the  joint 
bearing  and  shaft  flexibilities  (in  series) . 


The  transition  to  instability  occurs  as  the  rotor  speed  a>  *  .  As  o> 

passes  through  this  value,  one  of  the  roots  of  (38)  can  be  shown  to  cross  over 
into  the  half-plane  Re(A)  >  0,  the  crossing  taking  place  at 


K  m  io>Q  (41) 

corresponding  to  a  circular  whirl  of  angular  frequency  equal  to  the  critical 
frequency,  even  though  the  rotor  speed  is  twice  the  critical. 

For  co  >  2<oq  a  root  of  (38)  representing  an  expanding  whirl  results,  with 
the  frequency  of  whirl  increasing,  but.  slowly,  with  cd. 

If  is  allowed  to  approach  zero,  but  k  to  remain  finite,  then  cn^  in 
Equation  (41)  approaches  zero,  and  rotor  instability  at  all  speeds  results. 

This  result  agrees  with  earlier  studies  of  bearing-rotor  motion,  based  on  Equa¬ 
tions  (12,  (14),  Reference  (8). 

In  Equations  (28)- (41)  we  have  considered  motion  of  about  the  bearing 

center  0  for  a  journal  free  from  external  load.  Actually,  is  likely  to  vanish 

for  small  e,  judging  from  Figure  13,  and  some  external  load  Is  required  to  bring 

about  cavitation  on  the  expanding  film  side  and  generate  a  non-zero  F  and  k,  . 

r  b 

With  a  constant  external  load,  Equation  (33)  is  applicable,  with  £,  q  repre¬ 
senting  the  displacements  of  the  journal  center  from  the  equilibrium  portion 

(see  Figure  14) ,  and  F  ,  F  the  added  bearing  forces  due  to  this  added  displace- 

x  y 

ment.  At  the  same  time  e  must  be  sufficiently  small  to  justify  the  approxima¬ 
tions  in  Equation  (28)  to  Equations  (12)  and  (14). 

The  above  theory  historically  represented  the  first  successful  theoretical 
treatment  of  bearing  induced  whirl  (see  Reference  (j}  ),  and  it  was  obtained  at  a 
time  when  no  solutions  of  Equations  (9)  and  (19)  were  available.  By  increasing 
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with  e ,  one  obtains  a  higher  (dq  in  Equation  (40) ,  -thus  explaining  the  increased 
stable  speed  range  (39),  with  higher  static  loads.  Ye*.  a  more  sound  treatment 
should  be  based  on  solution  o£  Equation  (9)  with  cavitation. 

Such  solutions  have  been  carried  out  for  various  L/2a  (see  Reference  9  ) , 

so  that  the  functions  f  ,  f  and  their  derivations  at  6*  "0  are  available.  We 

r  t 

outline  the  small  displacement  stability  calculation,  based  on  them  and  on  Equa¬ 
tion  (27),  for  sizeable  £. 


On  Figure  14  ?et  0^  on  the  y-axis  be  the  equilibrium  position  of  the  journal 
center  for  a  given  static  load,  £,  T)  the  small  displacement  0, P  of  the  journal 
center  from  O^and  de,  da  the  corresponding  increment  in  e ,  a.  Equations  (27) 
for  small  de ,  da,  de ,  da  yield 

bf  bf 

dF  =  ip  a>  (-3—  de  +  de')  -  2daf 
r  de  de  r 

£f  bf  (42) 

dFt  =  in  cd  de  +  de')  -  2daft> 


where  and  and  their  derivatives  are  evaluated  at  e'  **  0.  In  rectangular 
components,  Equation  (42)  yields  for  the  added  forces  exerted  on  both  journals. 
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Attention  must  be  called  to  the  fact  that  dF^ ,  Y  have  opposite  sign  conventions, 
and  to  the  use  of  the  relations, 

da  51  £/Cae,  da  «  £/Cae  (44) 

in  obtaining  (43).  Equations  (44)  definitely  excludes  the  small  e  case. 

We  now  turn  to  the  solution  of  the  dynamical  rotor  equations  _ 

=  X,  mq  «  Y,  (45) 

considering  first  a  rigid  rotor  shaft.  We  assume  a  solution  ui  Equations  (43), 
(45)  varying  with  time  as  e^*" ,  and  put 


A.  =  VCD  ,  CD  =*  SXto 
o  o 


where  a>Q  has  the  dimensions  of  co, (1/sec)  ,  but  as  yet  undetermined ,  so  that  v,  s  are 
dimensionless.  There  results 


where 


A(v)  ~ 


sf  +  2 vf  +  Ee 
r  t 


sf.  +  2f  v 
t  r 


2  3  2 

men  v  mC  m  v 
o  o 


df  df 

_ r  _ t 

s  “57  +  v  “5? 

df  c3f 

S‘lf  +  V”dT’+  E 


0  (47) 


Equation  (47)  can  be  shown  to  be  stable  for  low  s  and  unstable  for  high 
enough  s.  The  transition  to  instability  occurs  by  V  crossing  the  pure  imaginary 
axis,  at  v  =  +iVQ  #  0  .  By  substituting  v  =  iv^  in  Equation  (47)  and  equating  to 
zero  the  real  and  imaginary  parts  separately,  two  equations  result.  The  second 


aquation  yields 


ft(0 


df  df  df  df 

(f  —4  ♦  £,  — r)  ♦  2  < f  —4  -  £  —4) 

r  b(  t  dc  t  dc  r  dt" 


>f 

_ EL, 

d«  '  ‘ 


(<<9; 


tnd  applying  this  to  the  first  ona,  one  oDtains 

of , 


-V  •  f2(£) 

2s 


(fr  '  efl(°]  “d7 


-  f. 


rdf 

cT  '  fi(£) 


s  Hi 
t  df ' 

of 

.  _ r 

'r  de ' 


(50) 


Dividing  Iquatlcn  (49)  by  (50)  and  recalling  (46),  thare  results 
t' 


_„3  2 

mC  cd 
o 


4f, 


n  La 


(51) 


Now  choose  o)q  so  that  e  *  1.  Then  Equations  (46)  and  (51)  lead  to 

2  2  4^  Hill 

“  ‘  mC3  f2(£) 


(52; 


while  Equations  (46)  and  (50)  yield 

K2  -  -  2 f 2 (e )  <do2.  (53) 

Equation  (52)  determines  the  critical  speed,  U)  ■  u>0  ,  at  which  rotor  motion 
becomes  unstable,  while  Bquation  (53)  yields  the  frequency  of  the  insipient  whirl 
at  that  speed  . 

The  above  theory  can  be  extended  to  include  a  finite  shaft  stiffness  coeffl- 
ciency  k,  by  properly  changing  the  definition  of  K.  (see  Reference  9  ) 

Turning  to  a  brief  consideration  of  compressible  lubricants,  it  is  evident 
that  cavitation  does  not  occur  since  a  gas  can  expand  to  fill  the  divergent  film 
half,  k  (p,  V)  curve  of  the  form  schematically  indicated  in  Figure  15  result*, 
where  V  is  the  specific  volume.  It  is  evident  that  equal  and  oppoalte  volume 
changes  +  £,V,  lead  to  unequal  pressure  changes  C5p)  +  ,  (£p)  .  Thus  mors  praaaur* 
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increase  is  obtained  on  the  convergent  film  half  on  Figure  12,  than  pressure  de¬ 
crease  on  the  divergent  film  half,  and  Equation  (5)  leads  to  a  non- zero  F  . 


The  Reynolds  Equation  (9)  is  replaced  by 

k  (h3c  It) +  is?  Ih3p  S? I  ■ 6,1  [-“  +  2bW  +  2h  It]  _  <54> 

Assuming  an  isothermal  (p,p)  curve  simplifies  Equation  (54)  to 

h  <h3p  ‘‘"Sae  <h3p  Tde)  "6p['aai^El  +  2pW  +  2h  a?3  (55) 

This  is  non-linear  in  p.  It  no  longer  involves  d,  a>  in  the  simple  form  a>  -  2d, 
as  in  Equations  (9)  and  (12)  .  For  stability  calculations  near  equilibrium 
positions,  solutions  varying  with  time  as  e^t,  require  further  solution  of  Equa¬ 
tion  (55),  due  to  the  term  dp/dt  on  the  right. 
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EFFECT  OP  DIFFERENT  BEARING  SUPPORT  STIFFNESS  CONSTANTS 


Thus  far,  the  bearing  supports  have  been  assumed  to  be  rigid.  We  now 
coraider  the  effect  of  resilience  or  flexibility  of  thase  supports  on  the 
ftablllty  of  the  rotors  and  shafts  considered  in  Sections  VI  -  VIII.  The  inertia 
of  the  supports  will  be  neglected. 

The  bearing  supportc  generally  posses  different  stiffness  constants  In 
two  mutually  perpendicular,  fixed  directions,  say,  the  horizontal  and  vertical 
directions.  Suppose  that  as  in  Section  VIII,  the  shaft,  too,  has  different  bending 
stiffness  constants  in  two  directions;  these  rotate  with  the  shaft.  The  net 
resilience  of  the  shaft  and  bearing  supports  turn  out  to  vary  with  the 
rotation  angle 

t  -  o)t  (1) 


whether  it  4s  exprt  f,ed  along  fixed  o  utating  axes.  One  is  thus  led  to 
linear  dilCarer*-'  -.’arions  or  rotor  motion  with  periodic  coefficients.  For 
the  roto:  u  ...,tlon  Vi  an  8th  order  system  with  periodic  coefficient  reeulte;  for 
special  cases  this  separates  into  two  4th  order  systems,  or  reduces  to  e  single 
4th  order  system.  For  the  rotor  of  Section  VI,  the  tranla.tional  motion  leads  to 
the  eyetem 


2  2 

px  +  (1+ofp  coe  2t)  +  (p  ein  2r)  »  0, 

dT  dT 


2  2 

d  6 

py  +  p  sin  2t  — y  4-  (1-o-p  cos  2t)  — j 

dr  dT 


(2) 


where  p,p,<J  are  proper  constants.  The  rotational  displacements  (with  unequal 
I,J,  leads  to  another  4th  order  system.  We  consider  the  stability  of  solutions 
of  (2). 

As  stated  in  Section  I  there  exists  special  solutions  of  linear  differential 
equations  with  periodic  coefficients  such  as  (2)  of  the  form 


_ .  .  At 

P(t )  e 


(3) 
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where  P  is  periodic  of  period  At  «  r.  and  X  is  a  constant,  four  such  solutions 
exist  for  the  4th  order  system  (2),  each  with  its  own  exponential  e  .  and 
general  solution  is  a  linear  combination  of  th'se  four.  For  stability  of  the 
motion  all  four  X's  must  be  in  the  half-plane 

Re  (X)  i  0  » 


If  any  one  of  the  X's  fails  to  satisfy  (4),  then  the  general  solution  of  Eq.  (2) 
will  become  infinite  at  t  increases. 


The  stability  conditions  (4)  are  similar  to  those  for  linear  differen¬ 
tial  equations  with  constant  coefficients;  the  difference  lies  in  the  great  dif¬ 
ficulty  of  obtaining  the  values  of  X. 


Since 

X(T+n)  Xt  Xu 

e  **  Ce  ,  C  “  e  , 


r  r.  \ 
/ 


the  solutions  (3)  have  the  property  of  being  multiplied  by  a  constant  C  when 
increases  by  the  period n  .  The  requirement  (4)  leads  to  the  condition  that  all 
four  C's  lie  within  or  on  the  unit  circle 


(6) 


The  derivation  of  (2)  proceeds  as  follows.  Denote  by  i  ,  1  unit  vectors 
along  the  (fixed)  (x,y)-axes  of  Fig.  1,  and  by  i  , ^  unit  vectors  rotating  with 
the  3haft,  and  corresponding  to  the  directions  of  the  principal  moments  of  Inertia 
of  its  cross  section.  A  force 


F  =«  F  i  +  F  J  **F(i  cost  -  j  ,?in  r)  F  (i  sin  t  +  j  cos  t)  (7) 

.xoyipx,  y 

applied  to  the  rotor,  will  produce  a  deflection  due  to  shaft  bending 

f1  (F  cos  t  +  F  sin  t)  i  +  f„  (-F  sin  T  +  F  cos  t)J  (8) 

i  x  y  c  x  y  ^ 

and  a  deflection 


f_  F  1  f  F  j 
3xo  4  y  o 


(9) 
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due  to  bearing  support  deflections,  where  f^,  are  the  flexibilities  (i.e.  re¬ 
ciprocals  of  stiffnesses)  of  the  shaft  in  the  i,  j  -  directions;  f  £  the 
flexibilities  of  the  bearing  supports  in  the  i  ,  j  directions.  Me  assume  that 

G  ^  O 

f^  >  f^,  f^^f^.  Adding  the  displacements  (8)  (9),  expressing  i,  j  in  (7)  in 
terms  of  i  ,  j  ,  and  introducing  the  notation 

-to  j? 


£»  -  <£1  +  £2  +  f3  +  V/2>  <■'  (£1-V/2V  0  -  (V£4>/2£n, 


(10) 


where  f  is  the  mean,  joint  flexibility  of  the  shaft  and  bearing  supports,  there 
results  for  the  net  rotor  center  deflection  x,  y 


x/f  *  (1  +  a  +  p  cos  2t)  F  +(psin  2  t)F  , 
m  x  y 


y/f^  *  (sin  2 r )Fx  +  (1  -  C  -  p cos  2t)  F  . 


(11) 


If  the  force  F  ,  F  is  due  to  the  rotor  inertia,  then 
x  y 

F  *  -m  ,  F  “  -m  . 


x  2  ’  v 

dL  y 


.  2 
at 


(17) 


Then,  replacing  F  ,  F  in  (11)  from  (12)  introducing  t  as  in  (1),  and  putting 
x  y  ^ 

2 


03 


2  1 


03 


o  mf 


■>  P 


(13) 


m 


03 


one  arrives  at  Eqs.  (2). 

Before  considering  the  solution  of  (2)  proper,  we  first  note  the  special 
cases  p  =  0,  and  0  =  0  of  equal  shaft  or  bearing  support  flexibilities.  In 
the  former  case 


p  "  °>  fr  =  f2  =  f8’ 


(14) 


the  trigonometric  terms  in  (2)  disappear,  the  periodic  terms  P(t)  in  (3)  reduce 
to  constants,  the  X' s  turn  out  to  be  pure  imaginary  and  correspond  to  sinusoidal 
solutions  in  x,  y  of  frequencies  03^ ,  03^  respectively  where 


>«i(i  2  +  13)  ’ 


03,, 


m(fg  +  f4) 


(13) 
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In  the  latter  case, 


°  -  °»  f3  “  f4  "  V 

by  using  rotating  axe6,  eqs.  (2)  may  now  be  reduced  to  Eqs.  6(7),  but  with 
k^,  k^  replaced  by 


■  +  fu> 


=  f  +  f, . 


k.  1  b’  k0  2  b 


Eqs.  6(13)  yields  the  unstable  range. 

"‘(Eb1+  f3)  “  n(thl+  V  ’*  (18: 

2 

while  for  oa  outside  this  range  the  motion  is  stable.  I:  erras  of  a,  P  Eqs.  (15) 

-  (18)  become 

P-  0;  u  *  (i  -  a),  u  ■  1  +  0,  (18] 

0  0;  1  -  P  <  n  <  1  +  P  .  (19; 

Summarizing,  it  is  seen  that  for  p  =*  0  the  constant  bearing  flexibility  may 
be  added  to  each  shaft  flexibility. 

For  general  a,  p  periodic  coefficients  persist  in  Eqs.  (2)  whether  fixed 
or  rotating  axes  are  used. 

Following  Floquet,  an  algebraic  equation  for  the  four  characteristic  con¬ 
stants  C  (see  Eqs.  (3)  -  (6)) can  be  formulated  if  four  (linearly  independent) 
solutions  of  (2)  are  available.  It  is  given  by  the  fourth-order  algebraic 
equation 

|Uij(vt)  -  C  5  |  =0,  (2o; 

where  j|  b^jj  is  fourth-order,  unit  matrix  (&„  =  1  for  i  =  j,  6^  =  0  for 

i  ^  j),  and 

I 

u  ii  =  xi('T)’  ui2  “  xi  ui3  =  ut4  “  yi'(T^  :  1  =  i-2-3*4  (21) 


are  four  solutions  of  (2)  whose  initial  values  are  equal  to  the  i-th  row  of  the 
unit  matrix  b„ .  however,  since  separate  integrations  of  (2)  over  the  interval 
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0  <  r  <  j;  have  to  be  carried  out  for  each  a,  p,  procedures  based  on  (20)  sppear 
economically  prohibitive,  even  with  the  aid  of  high  speed  computing  machines. 

Expansion  of  the  solutions  of  (2)  in  powers  of  P  is  a  possible  method  of 
integration  of  Eqs.  (2),  provided  both  A  and  P(t)  in  (3)  ai  e  expanded  in  powers 
of  p.  The  convergence  for  sizable  P is  very  slow,  however. 

Also  slow  is  the  convergence  of  the  "classical"  method,  due  to  A.  V.  Hill, 
babed  on  the  expansion  of  P(t)  in  (3)  in  Fourier  series  in  t.  A  significant 
improvement  in  its  convergence  was  effected  by  Foote,  Poritsky,  and  Slade,  in 
Ref.  (ID);  we  proceed  with  a  qualitative  description  of  this  method.  Ref.  11) 

First,  we  simplify  Eqs.  (2)  by  introducing,  as  in  Sections  III-V,  VII,  VIII 


Z  *  x  +  iy,  Z  -  x  -  iy,  i  «l/^T. 
One  obtains  from  Eqs.  (2) 

• •  2 ir  — 

Z  +  p  Z  +  (c  +  e  )  Z  =0, 

-  -  -2ir  " 

Z  +  P  Z  +  (o  +  e  )  Z  *  0. 


Setting 


_  AJ  niT  —  At  _  _  nix  .  .  .  _ 

Z=e  Z  A  e  ,  Z  *=  e  £  B  e  ;  n  ■  +  l,+3,. . 

n  n  -  — 


(22) 


(23) 


(24) 


substituting  (24)  into  (2)  and  equating  net  coefficients  of  each  exponential  to 
zero,  one  arrives  at  the  equations 


P  .  +  1  + 


(A+ni)' 


]  A 


+  a  B  =0 
n 


^  +  1  +  ..  .  2  B  +  p  A  _ 

n  (A  +  m)  J  n  n+2 


for  the  coefficients  A^,  B^,  where  n  =  +  1,  +  3,  ...  By  equating  determinants 
of  central  blocks  to  zero,  one  obtains 


A  (A)  =  0. 
n 


(26) 


The  limit  of  the  root::  of  (26)  as  n  gets  large,  yield  the  values  of  A. 


First,  it  will  benoted  that  renaming  the  (dummy)  summation  index  n  in  (24)  as  n  +  2 
or  n  -  2,  is  equivalent  to  replacing  k  by  \-2i,  A.+2i  respectively.  Hence,  tc  each 
value  of  k  .n  (3),  the  periodic  array  \+2ni  also  furnishes  possible  values;  they  all 
lead  to  the  same  value  of  C  in  (5)  since  e2ni=  1.  Next  note  that  changing  t  into 
-t  and  interchanging  Z,  Z  leaves  (2)  invariant.  Hence  the  four  A.'  s  are  arranged  in 
two  pairs,  +\^,  ‘  Finally,  from  (25)  it  will  also  be  seen  that  k  =  +  ni;  n=l,3... 

are  fourth  order  poles  Of  Ajj 00*  Now  a  function  of  X.  having  all  the  above  properties 
is  given  by 

,  .  .2  kn  .  .2  *1*  ,  ,2  lit  .  ,2  V\ 

(  stnh  —  +  sinh  —=*-)(  sinh  ~r~  -  sinh  — r— ) 

f(x)  =  - i - L- - f - (27) 

cosh  (n  X./2) 


It  is  therefore  to  be  expected  that  the  ratio  of  by  a  proper  function  of  n  ap¬ 
proaches  f(\)  as  n  — m. 


The  function  of f(\)  can  be  expressed  as  a  trigonometric  function  of 

(air/2)  ,  \  =  ai, 

and  its  numerator  as  a  quadratic 

x2  -  xOy-K.,  +  1)  +  K  =  0; 

where 

x  =  sinh2  =  (C-1)2/4C 

2  ^ln  2  *"2* 

=  sinh  — —  sinh  -y— 

2  ^ln  2  ^2rt 

K0  =  cosh  — y  cosh  '  -y—  . 

Thus  the  C’s  and  k'  s  can  be  found  once  ly  have  been  obtained. 


(28) 


(29) 


(30) 


The  convergence  cf  (the  roots  of)  A^K)  to  (the  roots  of)  f(\)  is  very  slow.  This 
can  be  demonstrated  from  the  convergence  of  the  infinite  products  in 


00  2 

cosf=  7T(1  - 

L  k 

k=l,3,... 


sm 


nx 

2 


OO 


TTu  -  M 

k=2,4,..k 


(31) 


to  their  limits.  In  Reference  (10  ),  K  ,  ty  are  computed  as  limits  of  certain 

sequences  of  K,  to  K, ,  and  K  to  K„  is  speeded  up  by  means  of  “convergence  factors" 
In  1  2n  / 
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which  are  basically  the  ratio  of  the  left-hand  members  of  (31)  to  their  finite  pro¬ 
duct.  sums.  For  details  the  reader  is  referred  to  the  reference  proper. 


Figure  lb  shows  the  unstac'-*--  regions  in  the  (p,p)  plane;  for  o  =  0,  to  the  right 
of  the  broken  lines;  for  o  =  .1,  to  the  right  of  the  solid  lines  with  a  triple 
zig-zag.  It  will  be  noted  that  the  single  unstable  range  for  0  =  0  breaks  up 
into  several  separate  ranges  for  o  =  0.1, for  small  p,  but  continue  as  a  single 
wider  range  for  larger  p.  For  large  0  ,  the  unstable  speed  ranges  separate  even 
further.  The  boundaries  of  the  unstable  region  correspond  to  passage  of  one  \-pair 
from  pure  imaginary  to  real  or  complex  values,  and  belong  to  the  loci 


\  -  0, 


-  -1, 


2  2 
X1  =X2 


corresponding  to  the  roots  of  Equation  29 
in  Figure  17. 

X1  =  °>  xi  -1-  x!  "  x2 


The  geometry  of  the  change  is  indicated 


i 

-I- 

T 


Fig.  17a 


Fig.  17b 


Fig.  17c 


\  2" 

On  Figure  16,  the  upper  and  lower  curves  correspond  Yo  X...  =0  and  Figure  17a 

\  ^  o  o 

and  "grow”  out  of  the  loci  (18)  for  p  -  0;  the  middle  one  Corresponds  to 
and. Figure  17c.  \ 


For  small  p.o  the  above  unstable  ranges  occur  near  p  =  1.  It  can  be  shown  that 
similar,  but  smaller,  unstable  ranges  occur  near 

2  2  2 
p  -  2  ,  3  ,  4 

that  is,  when  the  rotor  speed  is  near  a  submultiple  of  cno  7 


u)/o>  =  1/2,  1/3, 
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9 


V 


-  -)  T  I'  - . .  mriyiil  I  fcw  II  I  ,  »|  !04f  iff  li'ft  W„ll  1 


The  effect  of  static  friction  has  been  considered  in  Reference  (10);  it  narrows 
the  unstable  range  and  may  even  abolish  it. 

Similar  analysis  is  possible  when  unequal  bearing  support  flexibilities  are  added 
on  to  the  rotors  of  Section  VIII,  where  both  angular  and  tranalatory  motion  are  in¬ 
cluded.  In  particular,  the  (\,  ii)  motion  of  the  symmetrically  mounted  rotor  of 
Section  VI  considered  in  the  text  near  Equations  8(41),  8(42),  with  added  f^,  f ^ 
due  to  the  bearing  supports  (converted  to  an  equivalent  angular  stiffness  constant), 
leads  to  a  fourth  order  system  with  periodic  coefficients.  In  form,  these  equations 
are  the  same  as  for  a  rigid  shaft  passing  through  the  origin, with  different  angular 
stiffnesses  for  the  angular  deflections  X,  n.  This  case  has  been  studied  by 
Brosen's  and  Crandall,  Reference  (13,  using  the  same  method  of  solution. 

The  addition  of  f^,  f„  to  the  general  rotor  of  Section  8  leads  to  an  eighth  order 
system  with  periodic  coefficients. 

The  method  of  convergence  factors  is  applicable  to  higher  order  differential  equa¬ 
tions  and  will  be  discussed  in  another  paper. 
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XI 


APfMEP  K>MW»  QWm,  MMMKfiE. 

Thus  far,  the  rotors  have  been  assumed  balanced,  and  free  from  external 
forces.  We  now  consider  briefly  rotors  subject  to  gravity, and  unbalance  forces 
and  torques . 

With  the  addition  of  external  forces  (and  torques)  the  linear  equations  of 
motion  become  non- homogeneous.  Their  solution  consists  of  the  sum  of  a  particular 
solution  and  a  general  solution  of  the  homogeneous  equations  described  in  Sections 
II -VIII. 

Particular  solutions  will  now  be  considered  corresponding  to  a  constant 
force,  such  as  gravity  force;  then  to  the  unbalance  of  the  rotor. 

Under  a  constant  applied  force  (X,  0,)  such  as  the  gravity  force  X *■  -mg, 
eqs.2(l)  become 

mx  +  kx  -  X  ,  +  ky  -  0  (1) 

A  particular  solution  is  given  by 

x  ■  X/k,  y  •  0  (2) 

corresponding  to  the  static  shaft  deflection  under  the  action  of  this  force.  As 
stated  above,  to  this  is  added  the  general  solution  of  the  homogeneous  system 
2(1).  Hence,  the  motion  can  be  described  as  the  oscillation 2 (2),  not  about 
x  -  0,  y  -  0,  but  about  the  point  (2). 

Likewise,  for  eqs,  2(6),  the  solution  consists  of  a  sum  of  (2)  and  the 
damped  oscillatory  solutions  of  2(6). 

For  the  rotors  of  Sections  III  -  V,  a  constant  applied  force  producs  a  deflec¬ 
tion  in  the  direction  of  the  force,  as  well  as  at  right  angles  to  it.  In  par¬ 
ticular,  under  rotary  friction,  the  addition  of  the  force  (X,  0)  to  the  right- 
hand  members  of  eqs.  4(1)  yields  the  constant  displacement  given  by 
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X, 

0, 


(3) 


kx  +  cuy  » 
ky  -  f1  ok  ■ 


leading  to 


x 


kX 


k2  +  f  2 


tfX. 


2  2  2 
VT  +  (o 


producing  an  angular  displacement  relative  to  the  direction  of  the  force,  given  by 


tan  a  ■ 


(5) 


The  above  was  utilised  by  A.  L.  Kimball  in  measuring  the  coefficient  of 
Internal  damping,  reference  (13) ,  by  loading  a  horizontal  rotating  shaft  vertically, 
and  measuring  its  horizontal  deflection.  It  had  been  expected  that  a  would  vary 
with  the  rotor  speed,  and  with  the  load;  it  turned  out,  however,  that  within 
elastic  range  of  stress,  a  was  quite  constant,  and  small.  Thus,  if  a  rotary  fric¬ 
tion  model  is  used  for  internal  damping,  then  f^  must  be  assumed  to  vary  inversely 
as  a). 


Assuming  the  existence  of  a  normal  force  F  as  in  Section  III,  Eqs.3(2),  but 
changing  to  -k^,  one  finds  that  a  constant  applied  force  (X,  0)  produces  for 
small  k^/k  a  displacement 

X  klX 

x  -  £,  y  ■  — »  a  -  k^/k.  (6) 

The  elastic  energy  of  the  bent  shaft,  due  to  the  force  X  is  given  by 

V  -  X2/2k  (7) 

Since  the  shaft  is  revolving,  the  elastic  strain  of  any  element  of  the  shaft  can 
be  shown  to  vary  sinusoidally  with  time.  The  existence  of  the  force  F^  can  be  ex¬ 
plained  by  the  assumption  that  the  stress- strain  behavior  of  the  shaft  material  is 
described  by  narrow,  elliptical  hysteresis  loops,  and  this  assumption  leads  to  the 
energy  loss 
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V  -  ?taV 


(8) 


jtk, 


per  cycle.  Thu*,  a  i*  *  measure  of  the  "solid  friction*'  or  internal  damping 
of  the  material. 


For  the  rotor  of  Section  10,  a  constant  applied  force  Fx  produces  a  deflec¬ 
tion  given  by  eqs.  10 (  ) 


x  ■  f  F  (l+c  +  p  co*  2r) ,  y  ■  -f  F  p  ain  x  2r  (9) 

iq  x  in  x 

As  r  varie*  over  half  a  rotation,  At  ■  Jt ,  this  describes  a  circle,  show 
in  Fig.  8,  whose  diameter  is  at 

x  -  f_  F  (1+ff+p),  y  -  0,  (10) 

IQ  X 

*  "  -  P).  y  -  0, 


and  whose  radius 


We  now  consider  the  effect  of  unbalance.  Assume  the  rotor  symmetrically 
mounted  between  two  bearings,  as  in  Fig.  1,  and  suppose  that  the  unbalance  is 
likewise  symmetrical,  so  that  0^,  the  c.m.  of  the  rotor,  lies  in  the  midplane 
normal  to  the  z-axis,  at  a  distance  e  from  0^,  the  shaft  center. 


Denote  by  (x,  y)  the  components  of  the  displacement  of  0^,  from  the  unbent 
shaft  center  portion  0.  The  position  of  the  c.m.  0 2,  when  the  shaft  is  rotating, 

is 
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X  +  €  COS  CDt, 


(ID 


y  +  €  sin  cot 

assuming  that  as  t  ■  0,  0j02  Is  parallel  to  the  x-axis.  The  dynamical  equations 
2(1)  become 

2  2 
m  (x  +  e  cos  cot)  -  -kx,  m  — x  (y  +  c  sin  cut)  ■  -ky  (12) 

dtz  dt 

and  introducing  <dq  as  in  2(2), 

2  2  2  2 
5  +  o>Q  x  ■  e  co  cos  cot,  y  +  coQ  y  ■  e  cjo  sin  cot.  (13) 

The  effect  of  the  unbalance  is  thus  the  same  as  of  an  applied  rotating 
force,  equal  to  the  centrifugal  force. 


A  particular  solution  of  eqs.  (13)  is 


6  CO  COS  COt 

2  2 
0)  -  CO 

o 


co  sin  cot 
2  2 


Thus,  the  shaft  center  point  0^  describes  a  circle  of  radius  r,  where 


2 

co 

2  2 
CJJ  -  CO 

o 


and  00^  is  in  the  same  direction  as  0^  0^  if 


2  2 
co  <  o>o  , 


but  in  opposite  directions  if 


co  >  <or 


Eqs.  (14)  imply  that,  for  speeds  below  the  critical,  the  rotor  rotates  with 
"the  heavy  side  out";  at  speeds  above  the  critical,  it  moves  with  the  "heavy 
aide  in".  See  Fig.  19. 


A  plot  of  r/c  vs .  co  is  shown  in  Fig.  20,  Page  69. 
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y 


y 


£U  <  03 
O 


Fig.  19  '  1 

Fig.  20 

where  the  ordinate  plotted  represents  the  right-hand  member  of  (15),  but  without 
the  absolute  value  sign. 

From  the  above,  it  is  evident  that  the  critical  speed  co  can  be  defined  as 

o 

the  speed  at  which  an  unbalance  causes  a  larg^  (theoretically  infinite)  ampli¬ 
tude  of  shaft  displacement.  The  angle  between  OOj  and  O^Og  changes  from  0  to 
180°  at  the  critical. 

As  a)  becomes  infinite,  the  ordinate  of  Fig.  19  approaches  -1.  This  means 
that  0^  approaches  0  and  the  rotor  rotates  about  its  c.m. ,  which  approaches  the 
z-axls . 


At  oj  ■  o>o,  a  particular  solution  of  (2)  is  given  by 


A  e  a)Qt 

x  -  - 2 -  8ln  » 


«  “o' 

— r —  cos  co  t 
l  o 


This  may  be  described  as  motion  in  a  spiral,  known  as  "arithmetic"  or 
"Archimedean"  spiral,  with  the  radius  r  increasing  linearly  with  time,  indefi¬ 
nitely. 

If  a  "static"  friction  force  (-fx°,  -fy°)  is  added  to  the  right-hand  member 
of  (12),  then,  introducing  Z  -  x  +  iy  as  in  3(  ),  the  dynamical  equation  becomes 


'i+-Z+co2Z 
m  o 


2  Iwt 

e  e 
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This  leads  to  the  particular  solution 


-  2  lot 

Z  m  a)  e 

€  2  ,  ifo)  ,  2 

-0)  +  -  +  co 

m  o 


whence,  taking  real  and  imaginary  parts, 

x  ■  r  cos  (cot  -  $) ,  y  ■  r  sin  (cot  -  4>) 
where  (see  Fig.  20) 


2.2  .  ,2  2  .  211/2  ’ 

-  co  )  +  f  co  /m 

o  J 


2  2 

tan  <D  «  fco/m  (co  -  co  ). 


Fig.  21  '  Fig.  22 

A  schematic  plot  of  <J>  vs.  co  is  shown  in  Fig.  22  ,  and  r/e  vs.  co.is  in  Fig.  23 
There  is  now  a  gradual  change  in  from  0  to  -180°,  as  co  increases. 


Imfc  _  loot 

for  the  ratio  of  the  current  I  ■  Ie  to  an  applied  voltage  £  ■  E«  In  an 
R,  L,  C  circuit.  The  typical  resonance  curve  is  obtained  by  plotting  the  ab¬ 
solute  value  of  (24)  vs.  to.  Essentially,  (20)  differs  from  (24)  in  having  an 
extta  factor  ho  in  the  numerator.  The  maximum  of  Z/e  occurs  at  u>  ■*  cu>  >  co  , 

XB  O 

while  for  (24)  it  bccurs  *£,  <o  *  co  .  However,  if  f  co/o  is  small,  as  is  generally 

o  o 

the  case,  co  -  co  is  also  small, 
m  o 

Thus,  unbalance  causes  a  rotor  whirl  or  vibration:  ,  whose  amplitude  is 
large  near  the  speed  co  m  eoQ.  This  is  the  reason  for  the  term  "critical**  lot 
this  speed,  which  was  initially  defined  in  2(2)  as  the  vibration  frequency  of 
the  rotor  mass  m  vs.  the  spring  stiffness  k. 

For  the  rotors  of  Sections  3-5,  a  particular  solution  due  to  unbalance 

is  still  given  by  a  circular  motion  of  0^,  described  by  equations  similar  to 

(18)-(23),  but  with  r/e  and  0  given  by  different  functions  of  co.  This  is  the 

case,  even  when  the  solution  ~f  the  homogeneous  system  is  unstable,  as  for  rotors 

with  a  normal  force  F  as  in  ~q.  3(1). 

~n 

Addition  of  unequal  bearing  support  flexibilities  to  the  above  rotors  and 
shafts  yields  two  crltlcals  co^,  o^.  Now  unbalance  leads  to  an  elliptic  path 
for  0,  and  to  large  amplitudes  at  each  critical.  The  phase  changes  of  0  for  each 
one  near  co^,  co2,  have  the  effect  of  producing  for  a>j<  co  <  a>2  a  retrograde  motion 
of  0^,  that  is,  opposite  in  direction  to  the  rotation  of  the  rotor. 

Likewise,  for  the  symmetrically  mounted  rotors  of  Section  10  with  different 
shaft  and  different  bearing  flexibilities,  unbalance  leads  to  a  periodic  solution 
for  the  path  of  0.  For  small  p,  the  path  is  nearly  elliptical. 

Thus  far,  only  "translational  unbalance"  has  been  considered,  involving 
$c,  y)but  not  (\,  4).  We  now  turn  to  the  rotors  of  Section  VIII,  under  General 
unbalance,  where  the  rotor  displacement  Involves  x,  y  and  4,  and  the  unbalance 
involves  both  c.m.  off  the  z-axis,  and  the  principal  axis  of  inertia  off  the 
z-axls. 

Adopting  the  rotating  axes  of  Section  VIII,  we  suppose  that  the  rotor,  c.g., 
for  the  unbent  shaft,  is  located  off  the  z-axis,  at 
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X  ■  (  , 

X  * 


y  -  €, 


z  ■  z 


We  suppose  further  that  a  line  through  the  c.m.  parallel  to  the  z-axis  is  not  a 
principal  axle  of  Inertia  of  the  (undeflected)  rotor,  through  its  c.m.  Put 
Eq.  8(10)  in  the  form 


CD  ■  CD 


i  1,  +  +  CD  k-  -  (X.  -  cqi)  1 ,  +  (|!i  +  oA)  j2  +  cnk_  (26) 

x2r  y2^2  *2~  ^  ~Z  vZ 


Then  eq.  8(12)  is  replaced  *■;’ 


*  ■  l2  T1”,  -  Xx  y  “y2  ■  ly  ,  \] 
*  h  lA  y  %  +  J  “,2-  «  %  ] 


+  k,  (  -I  oj  -I  cd+Kcd 
2  L  X  Z  x2  y  *  y2  I 


where  Ix  y  ...  are  the  products  of  inertia  about  the  x2y2  ..  axes 

If  we  assume  "slight"  rotary  unbalances,  then  the  products  of  inertia 

I  ,  I  are  small  In  comparison  with  the  remaining  moments  end  products 
¥  i  y  z 

of  inertia. 

Substituting  for  i^,  J2,  k2  from  8(9),  differentiating  with  respect  to 
time  and  recalling  8(3),  one  obtains  i,  J  -  equations  which  differ  from  8(15) 


(a)  in  the  Introduction  of  added  terms  linear  in  x,  y,  4  on  the 

left  side  of  8(15)  involving  I  I._  ,  I  ,  and 

xy  xz  yz 

(b)  in  acquiring  the  right-hand  force  and  torque  terms. 


2  2  T  2  _  2 
■“  €x*  «“>  V  -  Xyz  “  *  '  Xxz  “  • 

As  regards  the  terms  on  the  left  which  are  linear  in  \,  4  and 

Involve  1,1,  these  may  be  considered  small  of  the  second 
xz'  yz 
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order,  end  will  be  neglected.  On  the  other  hand,  the  terms  involving  I  are 

xy 

equivalent  to  the  existence  of  an  angle  7  between  Ihe  principal  axes  of  the  rotor 
and  of  the  shaft  cross-section.  In  absence  of  Fig.  28,  one  is  led  to  a  form  of 
equations  equivalent  to  those  described  at  the  end  of  Section  VIXI,  involving  the 
angle  7. 

As  regards  the  terms  (28) ,  the  first  two  terms  represent  the  unbalance  forces  cor 
responding  to  the  righL-haad  members  of  (13);  the  last  two,  the  unbalance  moments 
due  to  the  angle  between  the  z-axis  and  the  nearest  principal  axes  of  inertia. 


Particular  solutions  of  the  resulting  differential  equations  can  be  obtained 
by  assuming  that  x,  y,  n  are  constants,  and  are  given  by 


x 


Nx  (uj.v) 
A(<d,v) 


V 


Nu  (m,v) 

A(oj,v) 


v 


0 


(29) 


where  N^...,  «r«  obtained  from  A  by  replacing  a  proper  column  of  A  by  the 

unbalance  terms  (28) . 


Since  the  axes  of  Section  VIII  are  rotating  axes,  the  solutions  (29)  represent 

a  circular  whirl  of  the  rotor  with  the  rotational  speed  cd,  with  the  unbalance 

forces  and  torques  balanced  by  the  elastic  forces  and  moments.  The  amplitudes 

2 

(29)  become  infinite  as  o>  approaches  the  "critical  speeds"  given  by  the  co  roots 

2 

of  A(co,0).  As  pointed  out  In  Section  VIII,  these  critical  speeds  given  by  the  co  roots 
as  transition  points,  at  which  the  solutions  of  the  homogeneous  system  8(15) 
change  from  stable  to  unstable  solutions.  £ Properly  speaking,  A  above  represents 
not  A  in  Fig.  8(16),  but  obtained  from  the  modified  system  resulting  from  the  replace- 

1 

ment  of  M  in  Fig.  8(14)  by  Fig.  (27). 
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XII 


VAEIABL&  SPEED 

Thus  far,  we  have  always  assumed  that  the  rotational  speed  of  the  rotor,  o>, 
while  arbitrary,  is  constant,  without  inquiring  how  this  speed  is  maintained.  Gen¬ 
erally,  the  needed  power  for  maintaining  this  speed  is  supplied  by  an  electric  motor, 
a  turbine,  a  diesel^or  some  other  "prime  mover".  It  has  thus  been  assumed  that 
these  "prime  movers"  supply  whatever  torque  is  needed  to  maintain  the  assumed  con¬ 
stant  rotational  speed,  to,  even  when  the  rotor  is  unstable  and  both  kinetic  and  po¬ 
tential  energy  is  required  for  the  increasing  amplitude  of  x,  y,  p. 

In  view  of  the  large  amplitude  vibrations  due  to  unbalance  near  the  critical 
speed  (See  Equations  XI (13) ,  XI (19)),  no  rotor  is  ever  designed  to  run  at  or  near 
its  critical  speed  if  it  has  one  such  speedy  or  if  it  has  several  criticals^ near 
any  one  of  its  criticals.  If  the  operating  speed  of  the  rotor  is  above  its 
critical(s),  then  it  must  pass  through  the  critical(s) .  For  a  perfectly  balanced 
rotor,  as  in  Section  II  the  rotation^and  the  translational  motion,  x,y^do  not 
interfere  with  each  other  even  when  the  rotor  speed  is  variable,  and  the  passage 
through  the  critical(s)  is  uneventful.  On  the  other  hand,  for  an  unbalanced  r<-  tor 
whose  speed  is  changing,  the  translational  motion  is  affected  by  the  variable 
speed.  We  proceed  to  consider  this  motion  of  an  unbalanced  rotor  whose  speed  is 
uniformly  accelerating  from  rest,  so  that  its  angle  of  rotation  9  and  speed  0  are 
given  by 

9  -  at~ / 2 ,  9  -  at.  (1) 

Unless  a  is  very  small,  the  solutions  of  Section  XI  say  Equations  XI (14) ,  are  quite 
inadequate  for  describing  the  rotor  behavior. 

Of  special  interest  is  the  rotor  motion  near  t  ■  t  ,  the  time  at  which  the 
rotor  speed  0  *  at  is  equal  to  the  critical  speed  cdo: 

t  *  a>  / a.  (2) 

o  o 

For  simplicity  we  consider  the  same  rotor  as  in  the  first  part  of  Section  11, 
that  is  a  centrally  mounted  rotor,  with  the  c.m.  O2  at  a  distance  e  from  the  shaft  axis 
01  ,  and  with  the  principal  inertia  axis  through  0^  parallel  to  the  n-axis.  For 

such  a  rotor  the  translational  motion  x,y  of  the  points  can  be  considered 

apart  from  the  angular  motion  (\,p)  of  the  axis,  and  we  proceed  to  consider  the 
former . 
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Denote  by  x,y  the  position  of  0^,  the  theft  center  point  in  the  plane  contain¬ 
ing  Og,  the  c.m.  Then  the  position  of  0 ^ »  assuming  that  0^  0^  is  initially  parallel 
to  the  x-axis,  is  given  by 

2 

x'  ■  x  +  €  cosO,  y'  ■  y  +  c  sin©,  0  ■  at  / 2. 

In  place  of  Equations  11(12)  one  now  obtains  the  dynamical  equations  of 
motion: 

,  2 

X  +  <D  X  " 

O 

..  .  2 

y  +  caQ  y  - 


d2  ,  at2. 
— ,  (€  cos  -r-) , 

dt4  L 

2  2 
d  /  ,  at  N 

— ,  fe  sin  — r— )  . 

dt^ 


(4) 


For  definiteness  we  consider  a  particular  solution  of  Equations  (4)  subject  to  the 
initial  conditions  of  an  unbent  shaft,  and  starting  from  rest: 


x  *  0,  x  ■  0,  y**0,  y-0  for  t  -  0. 


(5) 


The  general  solution  is  obtained  by  adding  the  solution  2(3)  to  the  above  particu¬ 
lar  solution.  (The  effect  of  gravity  changes  the  initial  value  of  x.) 


In  terms  of  xj  yj 

x  +  to  x'  -  to 
o  o 

2 

y'  +  co  y  "<D 
J  o  o 


and 

2 

€ 


2 

e 


the  coordinates  ol  the  c.m., 
2 

cos (at  /2) 
sin(at2/2) 


Equations  (4)  become 


(6) 


while  the  boundary  conditions  (5)  become 


x'  ■  £ ,  y'  *0,  x'  *0,  y'*0; 


these  can  be  reduced  to  (5) by 
troducing  the  complex  variable 


the  solution  x1  *  e  cos  u>ot,  y'  ■  0.  In- 


Z  -  x1  +  iy’ 


(8) 
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IV 


one  obtain*  from  (6) ,  (5) 


Z  +  o»o2Z  -  ®Q2  €  eiQt  11 ,  Z  -  0,  Z  -  0  for  t  -  0.  (S) 

The  solution  of  (9),  or  more  generally,  of 

u  +  ®  2u  *  f  (t)  ,  u  ■  0 ,  u  ■  0  for  t  »  0  (10) 

can  be  carried  out  in  terms  of  its  ’’Green.' s  Function"  G(t,s),  as  follows 
t 

u(t)-  /qG( t,s)  f(s)  ds,  (11) 

where  G  is  defined  as  follows: 

/  G(t,s)  ■  0  for  t  <  s, 


b2G  2 

— r  +  oj  G 

bfZ  ° 


0  for  t  >  s 


(12) 


dG(t,s) 

at 


S“t 


Physically,  G  represents  the  solution  of  (10)  for  an  "impulse"  force  at  the  time 
t  »  s 


f(t)  -  6  (t-s)  (13) 

where  5  is  the  "Dirac  Function",  while  Equation  (11)  represents  the  response  to  any 
f(t)  es  a  superposition  of  responses  to  equally  spaced  impulses. 

The  solution  of  (12)  yields 

sin  at  (t-s) 

G  (t,s)  =*  G(t  -  s)  *  - - -  for  t  >  s. 

o 

The  dependence  of  G  on  (t-s)  is  characteristic  of  linear  differential  equations 
with  constant  coefficients. 
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Applying  (11),  (14)  to  (9),  we  obtain 
Z(t)  ■  ®oc  JQ  ein  o»o  (t-e)  e1®*  ds 


(15) 


By  cxpreaeing  the  eine  in  term*  of imaginary  exponentials,  and  completing  the  squares 
of  the  quadratic  exponents,  the  (15)  can  be  put  in  the  form 


A(t) 


V 

2i 


^iflJbto/2  f.i^t  ;“toel^2/2du  .  e"l%t  ^to  eiOu2/2du 

L  -t  t 


(16) 


The  integral  in  (16)  is  related  to  the  Fresnel  integrals  defined  by 


E(v)  -  C(v)  +  iS(v)  -  L  einP  ' 2  dp, 


C(v)  »  /Q  cos  (np  12)  dp. 


S(v)  -  /Q  sin  (np  12)  dp, 


There  results 


Z(0 


o>  € 

o  e  o 


21 


—  |1  [E(v-v0,  -  E(-v0)]  -  e-“^  [«<*,.>  -  *<«„>]} 


(where 


(18) 


v  -  —  t,  v 
n  o 


«  t 
n  o 


Tables  of  the  functions  C(v),S(v)may  be  found  in  Jahnke-Ende.  The  function 
E(v),  plotted  in  the  complex  CH-iS  plane,  yields  the  well-known  Cornu  Spiral  indi¬ 
cated  schematically  in  Figure  24  .  For  large  v  ,  the  following  asymptotic  ex¬ 
pression  holds  : 


E(v) 


1+i 

2 


t  e 
n 


inv2/2 


+  .... 


(19) 


,  ,  1  1 

where  ....  refer  to  terms  involving  — y, 

v  v 
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%  u 


/VO. 

I  i 

■Jrn  i  i  i 

O-  '  / 


/ 

I  /  X  *N  \ 

I  ;  t  \  } 
x  v _ /  / 

The  path  described  by  ( is  given  by 


Fig.  24 


x*  +  iy'  ■  Z(t)  +  e  cos  <DQt 


and  aside  from  the  angle  6  ■  at  /2,  it  depends  on  the  angles 


•o  ■  Vo'  •<*.)  *  “o  12  • 

vhere  6  is  the  rotation  of  the  shaft,  9q  the  rotation  angle  of  a  rotor  rotating  with 
constant  angular  velocity  a>Q.  In  particular,  at  the  time  tQ,  v  ■  v  and  Equation 
(18)  yields 


ZU„>  -  j|  j«  ie°/2  E(v0)  +  .-3190«  (e(2vo)  -  *<v0)]J  .  (: 

The  maximum  of  jzj  occurs  generally  for  t  >  tQ. 

A  solution  similar  to  (11)  can  be  obtained  for  a  rotor  possessing  static  damp¬ 
ing,  as  in  Equations  2(6),  provided  G  is  replaced  by 

a  -  K<t->  .  +  b  .  „ 

D 

where  is  given  by  11(9),  Likewise,  the  integration  can  be  expressed  in  terms  of 
the  integral  E  but  for  non-real  arguments.  The  integral  E  is  related  to  the  error- 
function  but  again  for  complex  arguments. 


The  assumption  of  Equation  (1),  even  for  a  mas«ive  prime  mover  exerting  a  con¬ 
stant  torque,  is  not  necessarily  valid,  to  i;h.-  variable  torque  that  the  shaft 
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must  spply  to  maintain  the  motion  of  the  unbalanced  rotor.  Indeed,  the  kinetic  and 
potential  energies  of  the  rotor  are  given  by 


*  -  f  [ft 


(x  +  e  cosoj  +  (y  +  e  sin9)j^+  j  6 


Ji 


■  j  ^  +  y^  +  2e®  (y  cos®  -  x  sin®)j  +|  0' 


K  X2 


(22) 


„  k  ,  2  ,  2* 

V  »  j  (x  +  y  ) 


(23) 


The  first  two  Lagrange  equations  agree  with  (4)  while  the  third  one  yields 

t—  (y  cos®  -  x  sin®)  +  (me^+K)  ®j+  me0^(y  sin®  +  x  cos®)  ■  Q^,  (24) 

where  Q0  is  the  external  torque  about  the  z-axis.  If  the  above  solution  for'x.y,® 
is  valid,  then  the  torque-Q^  must  act  on  the  shaft.  In  general,  the  torsional 
stiffness  coefficient  of  a  shaft  is  of  the  same  order  of  magnitude  as  its  bending 
stiffness.  Hence,  the  assumed  Equation  (1)  must  be  corrected. 


If  the  prime  mover  is  not  very  massive,  the  effect  of  its  moment  of  inertia 
and  of  its  torque  vs.  speed  characteristic  must  be  taken  into  account. 
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XIII 


cagqgBP  mro? 

Thus  far,  in  Sections II  to XII,  many  basic  features  affecting  rotor  motion  and 
stability  have  been  considered.  For  clarity,  the  effect  of  each  complicating  fea¬ 
ture  was  treated  by  itself,  and  under  simplest  possible  assumptions.  Thus  the 
effect  of  bearing  lubricant  was  considered  in  Section  IX only  for  a  synsne  trie  ally 
mounted  rotor  consisting  of  a  single  disk  rotor,  but  it  was  neglected  in  Sections 
II  to  VIII,  X  to  XII. 

In  practice,  rotors  may  have  more  complex  structures,  and  many  of  the  features 
considered  in  Sections  II  tc  Xllmay  be  present  simultaneously.  As  an  example,  a 
steam  turbine  rotor  consists  of  a  fairly  massive  shaft  whose  thickness  varies 
axially,  and  on  it  are  shrunk  many  disks  carrying  steam  buckets.  This  "prime 
mover"  may  be  coupled  to  a  massive  turbine  generator.  The  "active  length"  of  the 
latter  forms  a  constant  radius  cylinder,  carrying  imbedded  current  coils  in  slots, 
and  for  a  two-pole  machine,  may  lead  to  different  bending  stiffnesses  in  two_ direc¬ 
tions.  The  coupled,  two- span  rotor  is  mounted  on  several  bearings  and  may  have 
an  overhung,  small  rotor  (voltage  regulator). 

A  complex  rotor  with  distributed  inertias  and  a  variable  thickness  shaft  is 
treated  by  breaking  It  up  into  many  sections,  each  one  with  its  own  translatory 
and  rotary  inertias  m^,  1^.  The  large  number  of  degrees  of  freedom  gives  rise  to 
a  higher  order  system  of  differential  equations,  with  more  criticals.  Inclusion 
of  damping,  bearings,  pedestal  support  flexibility  further  complicates  the  system, 
very  often  to  such  an  extent  that  high  speed  computing  machines  are  required  for 
their  treatment. 

To  illustrate  some  of  the  features  arising  from  an  increase  of  degrees  of 
freedom,  we  consider  a  rotor  with  two  disks,  as  shown  schematically  in  Figure  25  . 


Fig.  25 
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We  assume  that  the  shaft  and  the  disks  are  axially  symmetric,  and  that  each  disk  is 
perfectly  balanced. 


To  set  up  the  equations  of  notion,  -*e  start  with  the  displacement  load  equa¬ 
tions  similar  to  those  of  Section  VII. 


X1  "  *llFxl  +  a12Tyl  +  *13Fx2  +  a14Ty2  * 


^1  "  *21Fxl  +  *22Ty2  +  *23Fx2  +  *24Ty2  * 


x2  "  *31Fxl  +  a32Ty2  +  a33Fx2  +  *34Ty2  » 


^2  “  a41Fxl  +  a42Ty2  +  a43Fx2  +  a44Ty2  * 


*ij  "  *J1 


(1) 


Here  x^,  are  the  linear  and  angular  deflections  in  the  x,  z  plane  of'the  shaft 
at  the  center  point  n^;  x2 ,  the  corresponding  deflections  of  P2>  the  center 
point  of  m2>  while  Fxl  and  Tyl  are  the  external  force  in  the  x-direction,  external 
torque  about  a  line  parallel  to  the  y-axis,  applied  at  Fx2,  Ty2  a  similar  force 
and  torque  applied  at  P2.  Thus,  if  only  a  force  F^  is  acting,  then  the  deflect¬ 
ions  produced  by  it  at  P^,  P2  are  given  by  the  terms  in  the  first  column  on  the 
right  of  Equation  Q  and  so  forth. 

Solution  of  equations  (I)  yields 

1*1  -  Auxx  +  A12nl  +  A13x2  +  A14h2 


Tyl  "  A21X1  +  A22^1  +  A23X2  +  A24^2 


Fx2  "  A31X1  +  A32M1  +  A33X2  +  A34^2 


Ty2  “  A41X1  +  A42M1  +  A43X2  +  A44U2 


(2) 


where  A.,  is  the  "stiffness  matrix".  Both  the  A  and  the  flexibility  matrix  a.., 
ij  LJ 

by  Maxwell's  reciprocity  theorem,  are  symmetric: 


3ij  "  aJi’ 


Aij  "  Aji 


(3) 


I 
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tor  an  axially  symmetric  rotor,  tha  (y,\)  aquations  are  similar  to  Equations 
1  and  2,  and  can  be  obtained  from  thea  by  carrying  out  the  following  substitutions: 


J 


1 


(A) 


Consider  first  <»  ■  0,  that  is,  a  non-rotating  rotor.  The  dynamic  equations  of 
vibration  in  the  (x,s)-plane  are  obtained  from  Equations  (2)  by  replacing  its  left- 
hand  members  by  the  inertia  forces  and  torques  of  the  disks: 


Pxl  '  -1*1  •  V  ■  Px2  ■  -"2*2 •  ^2  ■  ‘Vl 


(5) 


There  results  an  eighth-order  system  of  four  differential  equations.  The  vibrations 
in  the  (y,s)-plane  follow  under  the  substitutions  (A) 

Assuming  a  solution  of  either  set  of  equations  in  which  time  enters  as  a  factor, 


(6) 


one  is  led  to  the  following  equation  for  p: 


yl 

X1 

y2 

X2 

^2 

-A2 

\rmip2 

A12 

A13 

A14 

A21 

A22-11P 

A23 

A24 

A31 

A32 

.  2 
A33'm2p 

A34 

A41 

A42 

A43 

aaa-v2 

where  the  first  row  on  the  top  Indicates  the  variable  occurring  in  each  column. 


(7) 


-82- 


The  (y,n)  frequency  equation  is  exactly  the  same,  provided  the  variables  be 
chosen  in  accordance  with  the  substitution  in  Equation  4,  hence  as  indicated  by 

the  second  row  above  the  determinant. 

Consider  next  id  >  0,  that  is  a  rotating  rotor.  For  it,  the  rotary  inertia 
terms  must  be  corrected  by  the  addition  of  gyroscopic  terms.  Recalling  Equations 
7(11)  and  7(12),  one  obtains  for  M^,  i L,  the  rate  of  change  of  moment  of  momentum 
of  each  disk,  ~  ** 


ji  -  (Il\  +  Vl*>  i  +  <¥l  '  Kl^l<0)  i 

*  (l2'^2  +  K2u2o>)  1  ■+  (I2u2  -  K2*2u>)  j , 


(8) 


and  the  i^- components  of  are  equated  to  T^,  T^2;  the  J-components  to 

to  -  ^xi'"Tx2‘  ^iere  results  now  a  system  of  eight  differential  equations  of  net 
order  16,  involving  all  eight  variables  x^,  . U2* 


It  is  possible  to  halve  the  order  of  the  system  by  introducing  the  variables 


Z^  *  x^  +  iy^,  -  i\^,  Z2  -  x2  +  iy2,  ^2  “  ^2  ”  ^2‘ 


(9) 


By  proper  manipulation  of  the  (x^ . M.2)-equations,  there  reaulta  a  system  similar  to 

the  (x,u)-equations,  except  for  the  addition,  in  the  principal  diagonal,  of  the 
terms 


0,  -kK^,  Oj-U^o*^. 


Under  the  assumption  :in  Equation (6) , these  added  terms  become 


(10) 


0,  K1U3pM1,  0, 


and  lead  to  the  net  inertia  terms 


(11) 


-m1p2Z_i,  (-I^2  +  K^cop)  M^,  -  m2p2Z2,  (-l2p2  +  K2ffip) 

There  result  in  the  variables  of  Equatiorg  (9)  f our  simultaneous  equations  which  ore 
the  same  as  the  (y  ,u) -equations .  The  determinant  of  the  coefficients,  when  equated 
to  zero,  resembles  Equation  7  for  a  non-rotating  totor,  provided  I ^ are  replaced  by 
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In  general,  for  only  a,  the  p-roots  are  real,  and  the  aolutlon  represents  a 
possible  circular  precession.  If  the  critical  speeds  are  now  defined,  as  in 
Section  12,  as  the  ay  values  for  which  the  precesslonal  frequency  is  equal  to  the 
rotational  speed  ®: 

P  -  ®;  (14) 

then  Equations  (13)  are  replaced  by 


* 


► 

* 


(15) 


Hence,  we  conclude  that  so  far  as  the  gyroscopic  effect  on  the  critical  speeds 
is  concerned,  the  static  vibrational  equations,  say  in  the  (x,z)-plane,  can  be  used, 
provided  the  rotary  inertias  of  each  disk  be  modified  as  in  Equation  (15). 

This  rule  applies  to  an  axially  symmetric  rotor  with  any  number  cf  disks. 

In  practice  even  the  calculation  of  the  stiffness  and  flexibility  matrices 
aiJ*  ^1J  l>ecome8  f^some  for  a  complex  rotor.  One  proceeds  by  step-wise  integra¬ 
tion  of  the  "beam  equations",  say  for  the  deflections  in  the  (x,z)-plane  for  x 
and  ii : 


dx 

dz 


El  ^  - 
dz 


dT 

“7^  +  S_  +  m 
dz  x  y 


0, 


(16) 


dS 


~r~  +  w 
dz  x 


for  the  deflection  x,  slope  u,  moment  T  ,  shear  S  ,  where  w  is  the  external  load 

y  x  x 


per  unit  z,  the  external  applied  moment  per  unit  z. 


Both  w,  and  m  are  derived 

f  y 


from  the  inertia  load  of  each  section,  and  are  concentrated,  say  by  being  divided 
equally  at  the  two  end  points. 

Assuming  a  value  of  cu,  one  solves  the  fourth-order  system,  subject  to  proper  boundary 
conditions  at  each  end.  For  a  single  rotor  span  of  length  £,  mounted  on  "fixed" 
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bearings  at  each  rod,  these  conditions  are 


x  -  0,  Ty  -  0  at  *  -  0, 

x  •  0,  Ty  -  0  at  z  -  t. 

On.  proceed,  by  applying  S,u.tlo»  (17)  and  ...unlng  eever.l  eet.  of  further  condition. 

the  vein.,  of  p  .nd  S£.  The  re.ultlng  .olutlon.  .re  then  interpol.ted 
to  eati.fy ,  ..,  x  -  0  et  .  -  t ,  .nd  Ty  .t  r  -  i,  1-  conputed.  Then  to  i.  varied 

until  T  vanl.hea  it  i  "  See  Reference.  (16) ,  (17) . 
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